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PREFACE 


A  comprehensive  program  has  been  undertaken  by  the  Federal  Aviation 
Administration  to  determine  the  suitability  of  the  Defense  Department's 
NAVSTAR/GPS  system  for  the  General  Aviation  (G/A)  user.  This  program  has 
in&luded  cost  studies  of  conventional  receivers  pertinent  to  the  G/A  user,  as 
well  as  design  studies  of  receivers  with  low-cost  potential  for  use  with  the 
standard  C/A  channel  satellite  transmission. 

As  part  of  this  overall  effort,  a  study  has  been  conducted  to  determine 
the  possible  simplification  in  receiver  complexity  for  the  G/A  user  resulting 
from  the  use  of  alternate  transmitted  waveforms.  Several  classes  of  alternate 
waveforms  were  considered  including  pulsed,  tone-ranging,  and  chirped-FM.  In 
addition,  the  present  waveform  with  substantially  increased  power  was  also 
considered  as  an  alternate. 

The  technical  analyses  carried  out  show  that  considerable  simplifications 
are  possible  with  both  approaches.  The  use  of  pulsed  waveforms,  combined  with 
emerging  matched-filter  technology,  makes  possible  a  greatly  simplified  receiver 
design.  If  increased  power  is  used  with  the  present  signal  design,  the  per¬ 
formance  will  also  be  substantially  improved;  alternatively,  performance  margin 
can  be  traded  for  simplified  receiver  Implementation. 
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SUMMARY 


*1.  Introduction 

Lincoln  Laboratory  was  tasked,  in  an  interagency  agreement  sponsored  by 
the  FAA  through  DOT/TSC,  to  examine  alternate  waveforms  to  the  presently 
planned  Clear  Acquisition  (C/A)  GPS  waveform.  The  purpose  was  to  review, 
select  and  propose  waveforms  which  would  perform  the  GPS  ranging  and  data 
carrying  functions  for  navigation,  using  the  GPS  satellites,  but  would  result 
in  greatly  simplified  (and  therefore  less  costly)  civil  user  receivers. 
Quantitative  performance  criteria  were  not  specified,  but  in  a  broad  compari¬ 
son,  performance  of  the  alternate  waveform,  in  terms  of  navigation  output, 
should  at  least  equal  the  performance  achievable  by  a  sequential  receiver 
using  the  GPS  C/A  waveform. 

Performance  criteria  to  be  included  were  ranging  performance,  data  hand¬ 
ling  capability,  time-to-first  fix,  acquisition  and  re-acquisition  capability, 
and  sensitivity  to  multipath. 

The  general  approach  specified  was  to  perform  a  quick  broadbrush  review 
of  waveforms,  including  the  present  C/A  waveform  with  higher  power  level,  and 
select  a  very  few  promising  waveforms  for  a  more  substantial  design  effort. 

The  more  detailed  designs  then  would  allow  credible  numerical  performance 
characteristics  to  be  obtained. 

Complexity  of  civil  user  equipment  was  to  be  the  major  criterion  for 
recommending  any  changes,  in  particular  alternate  or  additional  waveforms, 
to  GPS  because  user  equipment  cost,  which  reflects  the  complexity,  was  deemed 
to  represent  the  largest  component  of  the  total  cost  of  a  widely  accepted  and 
Implemented  system. 

The  principal  results  of  the  task  outlined  above  were  reported  in  brief¬ 
ings  to  the  sponsor  and  a  set  of  four  Working  Papers.  They  are  summarized  in 
this  report. 

2.  Outline  of  Principal  Results  and  Conclusions 

Following  a  short  review  of  the  principal  classes  of  waveforms,  based  on 
their  properties  to  determine  range  and  range  rate  when  optimally  processed, 
a  few  classes  were  identified  as  interesting:  the  burst  pulsed  high  peak  power, 
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low  duty  cycle,  waveforms  combined  with  some  method  of  bandspreading;  the 
pulsed  waveform  with  a  regular  PRF;  the  multitone  waveform;  the  present  C/A 
waveform  but  with  more  power.  For  each  class  a  detailed  design  was  proposed 
making  choices  for  a  number  of  additional  parameters  such  as  data  modulation, 
features  for  synchronization,  and  more.  While  there  was  often  a  good  reason 
for  specific  choices,  the  designs  and  their  characteristics  are  considered 
only  representative  of  what  can  be  obtained  using  the  waveform.  More  refined 
designs  may  result  in  different  and  more  pertinent  choices  of  parameters. 

Early  in  the  study  certain  "desirable"  signal  structure  features  were 
identified  which  were  adopted  in  the  designs  where  appropriate.  One  such 
feature  was  the  channel  sharing  method  based  on  TDMA  (Time  Division  Multiple 
Access).  The  pulse-burst  GPS  waveform  design  and  the  multitone  GPS  waveform 
both  use  TDMA. 

The  advantages  of  TDMA,  irrespective  of  the  specific  waveform,  are  rather 
compelling  for  a  low  cost  receiver  and  are  discussed  in  some  detail  in  Section 
7  of  this  report.  One  such  TDMA  format  would  allow  each  pair  of  antipodal 
satellites  to  transmit  in  one  of  twelve  timeslots,  the  total  cycle  of  twelve 
slots  with  fixed  assignments  repeating  itself  indefinitely.  Each  timeslot  can 
be  further  subdivided  into  signal  segments  performing  specific  functions  and 
into  a  "guardtime"  to  guarantee  that  signals  transmitted  in  consecutive  timeslots 
by  satellites  at  large  range  differences  to  a  user  can  not  overlap.  In  the 
pulse-burst  design  a  first  segment  consists  of  a  burst  of  pulses  for  synchroni¬ 
zation  and  ranging.  A  second  segment  of  pulses  transmitted  at  a  lower  duty 
cycle  serves  the  purpose  of  transmitting  the  data.  In  the  multitone  design 
the  first  segment  is  reserved  for  transmission  of  pure  carrier  for  phase-lock¬ 
ing,  a  short  second  segment  is  reserved  for  transmission  of  tones  for  ranging 
for  which  phase  locking  can  occur  very  fast  since  the  carrier  frequency  offset 
is  known  from  the  first  segment,  and  in  a  third  segment  DPSK  modulated  carrier 
is  radiated  for  data  transmission.  Briefly  stated  the  principal  conclusions 
on  the  pulse-burst  design  (Section  4  in  this  report)  are  as  follows. 


If  a  high  peak  power  (2kw) ,  low  duty  cycle  pulse  burst  waveform,  having 
about  one-third  the  average  power  and  the  same  information  content  as  the 
current  C/A  waveform,  were  added  to  or  would  replace  the  C/A  signal  on  the  GPS 
satellites,  it  would  be  possible  to  produce  a  high  performance  navigation  re¬ 
ceiver  of  considerably  less  complexity  than  the  present  single  channel  sequen¬ 
tial  receiver.  The  specific  PN-encoded  pulse  burst  waveform  design  and  the 
receiver  design  take  maximum  advantage  of  the  features  of  the  emerging  CCD 
(Charge  Coupled  Device)  matched  filter  technology.  The  receiver,  as  compared  to 
a  receiver  using  the  C/A  signal,  exhibits  comparable  or  better  accuracy,  greater 
update  rate  (of  the  order  of  one  per  second)  which  can  be  made  adaptive  depend¬ 
ing  on  user  platform  maneuver  condition,  and  lower  time  to  first  fix  (of  the 
order  of  20  seconds) .  This  receiver  would  exhibit  less  sensitivity  to  multipath 
and  to  signal  dropouts  because  it  would  continuously  track  all  satellites  in 
view  and  could  essentially  instantly  re-acquire  a  lost  signal. 

These  advantages  of  an  additional  waveform  must  be  weighted  against  the 
cost  and  difficulty  of  producing  high  power  pulses  on  board  the  satellites. 

A  brief  survey  has  indicated  that  suitable  space  qualified  power  generation 
hardware  is  or  soon  will  be  available.  Also  interference  and  interoperability, 
even  if  the  new  waveform  had  to  replace  the  C/A  waveform,  do  not  seem  to  present 
insurmountable  problems.  No  detailed  analysis  however  was  done  on  these  issues. 

The  advantages  of  an  additional  waveform  consisting  of  multitones  were  not 
as  clear-cut.  The  specific  design  (Section  5  )  allowed  demonstration  of  the 
capability  of  adequate  performance  with  considerably  less  peak  power  require¬ 
ments  (3  dB)  than  the  pulse  burst  system  described  above.  The  system  would  be 
sensitive  to  interference  however.  The  core  of  the  receiver  hardware  would 
consist  of  a  bank  of  PLL  (Phase  Locked  Loops)  operating  in  a  narrow  band  around 
a  set  of  known  but  changing  frequencies.  Unless  great  strides  are  made  in 
developing  cheap  PLL's  in  the  half  megahertz  range,  no  great  cost  savings  are 
expected. 
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A  third  class  of  waveforms  which  was  considered  (Section  6)  consists  of 
an  unending  PN-encoded  pulse  train,  with  a  PRF  (Pulse  Repetition  Frequency) 
specific  for  each  antipodal  satellite  pair.  The  advantage  of  this  waveform 
over  the  TDMA  pulsed  GPS  design  lies  in  the  fact  that  a  lower  duty  cycle 
satellite  transmitter  is  needed  and  the  user  clock  stability  requirements  can 
be  relaxed  (from  one  to  6  parts  in  10^  per  second).  As  was  the  case  for  the 
TDMA  designs,  the  use  of  signals  from  more  than  the  minimal  four  satellites  is 
possible  at  little  additional  computational  cost  while  eliminating  problems 
associated  with  satellite  selection  algorithms  and  problems  associated  with 
change-over  between  selected  configurations.  The  specific  PRFDMA  (Pulse  Re¬ 
petition  Frequency  Division  Multiple  Access)  design  requires,  during  a  signal 
acquisition,  additional  computation  for  pulse  sorting  which  was  absent  in  the 
other  designs.  And  finally,  the  class  of  the  present  CH  PN  encoded  GPS  wave¬ 
form  with  a  slight  parameter  variation  was  considered.  The  waveform  is  ill- 
suited  to  a  sequential  receiver  for  the  data  collection  function.  To  resolve 
this  difficulty  satisfactorily  one  may  resort  to  a  separate  data  collecting 
channel  and  an  additional  signal,  or  dedicate  time  intervals  to  collecting  data 
while  suspending  the  tracking  function.  In  general,  more  power  yields  better 
performance  of  the  delay  locked  loops  and  phase  locked  loops  resulting  in 
faster  (re-)  synchronization  and  possibly  higher  ranging  accuracy.  The  higher 
power  level  makes  matched  filter  applications  attractive.  The  CCD  correlator, 
around  which  the  burst-pulse  GPS  waveform  was  designed,  is  less  well  suited 
for  the  C/A  waveform  because  of  the  long  code  (1023  chips)  and  the  need  to  make 
it  programmable,  although  technology  advances  may  change  this  situation. 

In  this  study  we  also  investigated  several  important  issues  transcending 
specific  waveform  choices,  such  as  handover  capability  from  an  alternate  wave¬ 
form  to  the  P-mode  and  the  clock  stability  requirements  for  several  range  up¬ 
date  rates.  The  former  was  found  not  to  present  any  serious  difficulties.  The 
clock  stability  requirement  for  the  TDMA  pulse  (burst)  design  and  the  multi- 
tone  design,  apart  from  bias  and  frequency  offset,  was  characterized  by  (the 
square  root  of)  the  Allan  variance  of  one  part  in  10^  per  second  ,  to  keep  the 
clock  error  contribution  to  the  ranging  error  below  13  feet.  For  the  PRFDMA 

design  the  clock  stability  requirements  could  be  relaxed  to  about  6  parts  in 
7 

10  for  the  same  contribution  to  the  ranging  error. 
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1.  INTRODUCTION 


1 • 1  Background 

The  basic  reason  for  studying  the  possible  use  of  GPS  by  civil  aviation 
is  that  the  system  will  almost  certainly  become  fully  operational  and  the  poten¬ 
tial  for  taking  advantage  of  this  navigational  resource  by  non-military  users 
will  also  exist.  The  system  is  being  developed  and  will  be  maintained  by  DOD, 
under  present  plans;  hence  a  navigational  aid  with  outstanding  capabilities  in 
accuracy  and  coverage  will  be  in  place  for  use  by  the  military.  This  system 
can  also  be  used  by  the  civilian  sector,  as  OMEGA,  LORAN  and  TRANSIT  have  been  used. 

A  difficult  issue  to  address  is  the  establishment  of  cost- 
performance  goals  for  airborne  navigational  equipment.  It  is  not  unreasonable 
for  a  user  to  expect  a  new  generation  of  equipment  to  provide  better  performance 
at  a  lower  coat  than  the  one  it  replaces,  particularly  if  the  major  development 
costs  of  the  system  (including  some  classes  of  user  equipment)  are  borne  by 
the  Defense  Department.  Performance,  however,  is  a  multidimensional  quantity 
which  includes  not  only  position  determination  accuracy  and  update  rate,  but 
coverage,  reliability  and  the  efficiency  with  which  the  user's  interface  with 

the  equipment  tells  him  what  he  wants  to  know  regarding  where  he  is  and  how 
he  is  to  get  where  he  wishes  to  go.  Without  performance  criteria,  new  systems 
concepts  tend  to  be  compared  with  existing  systems  in  artificial  ways,  such  as 
how  much  better  the  new  one  will  be  with  respect  to  isolated  performance 
measures.  What  one  wants  to  know,  instead,  is  how  ouch  better  off  the  user 
will  be.  In  a  comprehensive  and  general  sense,  if  he  adopts  the  new  system. 

This  viewpoint  does  not  lend  itself  to  numerical  precision,  but  we  intend  to 
be  guided  by  it  in  assessing  the  various  alternatives  and  their  impact  on 
user  equipment  within  the  scope  of  this  study. 

1.2  Scooe  of  Study 

The  present  study  is  focused  on  user  equipment  suitable  for  civil  air¬ 
craft,  and  using  GPS  signals.  We  are  primarily  concerned  with  receiver  sim¬ 
plicity,  reliability  and  operational  characteristics.  No  attempt  will  be  made 
to  predict  actual  dollar  costs  of  equipment  alternatives,  since  predictions  of 
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this  kind  require  in  each  case  a  detailed  design  and  production  cost  breakdown* 
together  with  a  practiced  judgment  regarding  market  realities  which  we  do  not 
pretend  to  possess.  However,  we  are  very  much  aware  that  cost  to  the  user  in 
dollars  is  the  real  issue,  and  we  intend  to  address  that  issue  as  closely  as 
possible  by  using  engineering  judgment  in  place  of  the  techniques  of  production 
cost  forecasting.  We  also  believe  that  this  is  a  reasonable  approach  to  the 
problem  at  its  present  state,  since  the  system  alternatives  to  be  considered 
are  sufficiently  varied  as  to  cause  the  resulting  receiver/processors  to 
differ  fundamentally  from  one  another.  At  a  later  stage,  when  design  choices 
of  a  much  more  specific  kind  are  to  be  made,  within  the  context  of  a  chosen 
signal  and  data  format,  then  hard  dollar- cos ting  efforts  would  be  appropriate. 
But,  in  this  study  we  rely  heavily  on  complexity  as  an  indicator  of  cost,  along 
with  a  general  knowledge  of  the  cost  trends  of  the  major  component  classes 
which  are  bound  to  be  represented  in  the  user  equipment. 


There  is  need  for  a  study  of  equipment  alternatives  for  civil  aviation 
users  even  in  the  context  of  the  presently  planned  waveforms  and  signal  para¬ 
meters.  The  military  themselves  envision  an  entire  spectrum  of  user  equipments 
for  these  signals.  The  questions  of  the  relative  advantages  of  multichannel 
versus  sequential  receivers  and  the  possible  use  of  the  wideband  P  channel  ver¬ 
sus  exclusive  use  of  the  relatively  narrowband  C/A  channel  are  still  open,  in 
the  context  of  minimum-cost  receivers  meeting  civil  performance  goals.  More 
detailed  questions,  such  as  loop  parameter  optimization  and  processing  algo¬ 
rithms  for  position  determination  and  satellite  selection  also  merit  study. 

However,  the  problem  is  much  broader  than  this,  since  it  includes  the 
general  option  of  another  signal  for  civil  users.  No  detailed  ground  rules 
have  been  stated  regarding  such  a  signal,  beyond  a  general  willingness  to  pro¬ 
vide  a  modest  amount  of  space,  weight  and  power,  aboard  the  satellites,  for  this 
purpose.  Obviously,  such  an  add-on  package  must  not  interfere  with  the  original 
equipment  in  any  way.  Still  another  option  is  the  modification  of  one  or  both 
of  the  GPS  signals  to  facilitate  the  design  of  low-cost  receivers  for  civil  use. 
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It  Is,  therefore,  natural  to  divide  the  study  Into  two  parts:  one  devoted 
to  exploring  the  potential  of  parameter  variations  within  the  framework  of  the 
present  GPS  signal/data  structure,  and  the  other  devoted  to  an  evaluation  of 
the  possibilities  of  other  generic  classes  of  waveform  and  data  structure. 

The  problems  of  carrier  frequency  and  average  power  are  common  to  both  parts 
of  the  study,  while  consideration  of  short-term  and  long-term  duty  cycles  for 
the  satellite  transmitter  also  arise  in  the  second  part. 

The  signal  requirements  for  a  one-way  ranging  system  such  as  GPS  are 
basically  a  time-tagged  ranging  waveform  and  a  means  of  conveying  data  which 
allow  the  user  to  compute  the  satellite's  position.  The  time  tagging  consists 
of  the  provision  of  recognizable  points  in  the  waveform  together  with  an  in¬ 
clusion  in  the  data  of  the  times  of  transmission  of  each  of  these  points.  The 
present  GPS  waveforms  (P  and  C/A)  are  CW,  using  biphase  PN  coding  for  ranging, 
and  synchronous,  biphase  DPSK  modulation  for  data  transmission.  The  major 
alternative  ranging  waveforms,  especially  pulsed  waveforms  and  CW  tone  ranging 
waveforms  will  be  considered.  Each  ranging  waveform  can  be  combined  with  any 
of  several  feasible  data  modulation  techniques  to  complete  a  waveform  definition. 
In  addition,  the  data  structure  can  be  varied.  This  is  an  important  consider¬ 
ation  since  the  relevant  data  remains  unchanged  for  long  periods  in  satellite 
memory  and  can  be  transferred  to  the  user  for  acquisition  of  precise  and 
approximate  information  on  the  ephemerides  of  all  the  satellites,  clock 
anomalies,  ionospheric  propagation  models,  and  so  on. 

Besides  these  alternatives  for  the  signals  from  individual  satellites, 
there  are  important  alternatives  for  sharing  the  available  channel  among  the 
24  satellites.  Presently,  GPS  uses  code  division  (CDMA),  which  combines  well 
with  the  PN  ranging  waveform.  The  other  principle  alternatives,  TDMA  and 
FDMA,  will  be  considered  for  each  of  the  waveform  alternatives. 
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Thus,  the  signal  structures  to  be  studied  can  be  broken  down  in  a  number 
of  ways: 

GPS  modification  or  new,  additional  signal 

Carrier  frequency  and  power  levels 

Ranging  waveform  used 

Data  modulation  technique  employed 

Structure  of  the  transmitted  data 

Nature  of  the  multiple-access  technique  used. 

Each  of  the  many  alternatives  resulting  from  this  multidimensional  break¬ 
down  will  be  examined  from  the  point  of  view  of  the  implications  for  the  re¬ 
ceiver  in  terms  of  complexity,  as  discussed  above. 

1.3  Methodology 

It  would  be  desirable  to  establish  the  potential  of  each  of  the  signal 
options  described  above  for  yielding  satisfactory  performance  using  a  simple 
receiver.  As  mentioned  earlier,  this  performance  is  not  a  single  number,  but 
a  set  of  qualities,  many  not  amenable  to  precise  numerical  measure.  Some 
performance  features  which  can  be  assessed  numerically  are: 

Accuracy  of  a  single  position  fix 

Update  rate  of  effectively  Independent  fixes 

Information  data  rate 

Time  to  a  first  fix  (TTFF)  from  start  with 
varying  degrees  of  a  priori  information. 

Even  these  qualities  are  described  by  averages  and  probability  distributions 
because  they  depend  upon  many  variables,  such  as  satellite  geometry  and  user 
latitude.  Other  qualities  must  be  described  parametrically,  such  as: 

Probability  of  loss  of  ranging  data  and  time  to 
re-acquisition  as  a  function  of  user  dynamics 
and  orientation. 
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A  final  class  of  performance  features  is  qualitative,  but  no  less  important, 
in  nature.  These  include: 

The  initialization  requirements  of  the  system, 
i.e.,  initial  loading  of  current  GPS  system 
data  (if  this  is  required)  and  the  ease  with 
which  this  is  accomplished. 

The  nature  of  the  output  information  displayed 
to  the  user  and  the  ease  with  which  he  can 
obtain  navigational  guidance  from  his  equipment. 

Fortunately,  these  last  two  items  are  largely  independent  of  the  signal  wave¬ 
form  and,  although  somewhat  dependent  on  data  structure,  can  to  a  large  extent 
be  discussed  as  features  of  any  user  system  based  on  GPS,  as  compared  to  those 
based  on  VOR/DME,  Loran-C,  etc. 

Returning  to  the  problem  of  complexity/performance  goals,  these  are  of 
direct  relevance  in  the  case  of  two  of  the  performance  qualities  which  can  be 
most  reliably  determined  for  a  system  concept,  namely  positional  accuracy  and 
update  rate.  These  requirements  vary  greatly  between  the  air  traffic  control 
extremes  of  crans-oceanic  flight  and  final  terminal  approach.  Correspondingly, 
the  cost-effectiveness  of  a  navigational  system  depends  more  on  which  of  these 
ATC  regimes  it  serves  satisfactorily  than  whether  its  positional  accuracy  is 
100  feet  or  1000  feet.  To  relate  the  two  requires  the  missing  ingredient  of 
performance  goals. 

If  every  signal  option  were  evaluated  in  terms  of  all  the  performance 
measures  discussed  above,  a  large  matrix  of  numerical  and  qualitative  data 
would  be  obtained.  Each  element  of  that  matrix  would  be  the  result  of  optimiz¬ 
ing  a  receiver  design,  balancing  the  various  conflicting  performance  require¬ 
ments,  for  that  particular  signal  option.  Besides  being  a  task  of  impractical 
magnitude,  the  results  would  be  only  as  good  as  the  optimization  performed, 
and  the  subjective  choices  made  in  such  an  optimization. 


Instead  of  such  a  thorough  evaluation  we  divided  the  waveforms  up  Into 
major  categories  based  on  the  prime  function  of  the  GPS  waveform,  namely,  rang¬ 
ing  and  by  the  exercise  of  judgment  rather  than  detailed  evaluation  rejected 
those  which  did  not  seem  competitive.  This  process  led  to  three  major  classes 
of  waveforms,  namely  the  pulsed  TOMA  waveform,  the  tone  ranging  TOMA  waveform 
and  the  present  CW  class.  Rather  than  consider  all  possible  variations  and 
refinements  in  each  class,  a  complete  system  design  was  performed,  based  on 
a  large  set  of  specific  choices,  for  each  major  class.  The  resultant  candidate 
is  then  used  as  a  vehicle  for  performance  analysis  and  in  its  broad  conclusions 
representing  the  whole  class.  If  an  actual  system  design  choice  were  to  be 
made  one  could  further  investigate  the  variations  in  each  class  for  optimiza¬ 
tion.  Such  work  was  considered  beyond  the  scope  of  this  study. 

Some  issues  were  found  to  transcend  the  choice  of  waveform  modulation. 

For  example,  it  appears  that  multiple  access  by  time  division  is  inherently 
desirable,  from  the  point  of  view  of  receiver  costs,  because  one  receiver/ 
processor  channel  can  then  perform  ranging  on  all  visible  satellites  in  turn. 

A  less  obvious  result  is  that  with  ranging  data  available  from  all  usable 
satellites,  an  elaborate  satellite  selection  algorithm  can  be  dispensed 
with,  and  superior  navigation  obtained  with  a  nominal  increase  in  the  com¬ 
plexity  of  the  position  computation  algorithm.  Some  general  conclusions  such 
as  this  can  be  made  without  detailed  system  designs,  but  often  the  basic 
truths  of  a  subject  as  complex  as  this  one  are  perceived  only  after  working 
out  many  detailed  examples. 

Another  example  may  be  the  role  played  by  multipath.  It  is  hard  to 
conclude  much  from  general  statements  on  multipath,  but  an  assessment  of  its 
effects  on  a  few  candidate  signal/receiver  designs  may  permit  some  generali¬ 
zations.  It  was  decided  not  to  pursue  the  effects  of  multipath  in  our  detailed 
designs. 
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The  main  results  of  this  study  are  the  selection  of  a  few  promising 
alternate  waveforms  for  which  a  GPS  waveform  and  receiver  design  was  carried 
through  the  initial  stages  to  verify  performance  claims.  It  emerges  that  the 
class  of  pulsed  TDMA  waveforms  is  promising  enough  to  warrant  further  develop¬ 
ment  and,  perhaps,  detailed  further  work  by  the  civilian  agencies  in  order  to 
exploit  the  potential  of  NAVSTAR  GPS. 

Most  of  the  material  presented  in  this  work,  and  portions  not  reported 
here,  has  been  published  in  the  form  of  Lincoln  Laboratory  working  papers. 

1.4  Further  Outline  of  the  Report 

In  Section  2  we  discuss  the  winnowing  of  the  options  based  mainly  on 
the  capability  to  perform  ranging,  which  is  still  the  most  important  function 
to  be  performed.  Out  of  this  discussion  emerge  signal  design  concepts  qualify 
ing  for  further  evaluation.  They  are: 

a)  Section  3,  the  basic  GPS  C/A  signal  structure, 
but  enhanced  in  power  level  by  a  significant 
amount. 

b)  Section  4,  a  TDMA  pulsed  waveform  system 

c)  Section  5,  A  TDMA  tone  ranging  system. 

d)  Section  6,  a  PRFDMA  pulsed  waveform  system. 

In  Sections  7,8  and  9  we  discuss  issues  transcending  specific  waveform 
choices  such  as  the  operational  implications  of  TDMA,  handover  capability  and 
the  role  of  clock  stability  in  a  low-cost  receiver. 


2.  REVIEW  OF  WAVEFORM  CHOICES 


2.1  Choice  of  Waveforms  for  Ranging 

In  this  section  we  review  quickly  the  several  ranging  methods  and  provide 
some  rationale  for  having  limited  our  study  to  a  few  classes  of  waveforms  for 
each  of  which  we  will  design  a  prototype  GPS  system. 

For  purposes  of  navigation  the  user  Is  In  an  unknown  position  (within 
bounds)  and  tries  to  establish  his  position  relative  to  several  sources  at 
known  locations  (with  the  sources  often  providing  the  data  from  which  their 
location  can  be  determined) .  One  method  of  obtaining  relative  location  is 
through  ranging  or  range  differencing,  and  then  establishing  the  user  position 
by  multilateration. 

Some  navigation  systems,  such  as  IMS,  do  not  rely  on  external  signals, 
but  are  capable  of  integrating  displacement  and,  given  correct  initial  condi¬ 
tions,  provide  user  position.  We  do  not  consider  such  systems,  not  even 
as  ranging  aids ,  but  concentrate  instead  on  methods  of  ranging  based  on  exter¬ 
nally  transmitted  and  received  signals. 

2.2  Basic  Types  of  Ranging  Systems 

There  are  many  possible  system  approaches  to  ranging.  First  there  is  the 
interrogate-reply  type,  where  the  user  interrogates  several  sources  (e.g.  satel- 
lites)  and  deduces  two-way  range  from  the  delay  in  the  reply.  The  most  rudi¬ 
mentary  form  of  interrogate-reply  is  the  radar  system  based  on  an  interrogation 
with  RF  pulses.  If  an  object  reflects  a  portion  of  the  energy  in  the  pulses, 
the  delay  with  which  the  echo  is  received  is  proportional  to  the  two-way  range. 
Another  form  is  the  cooperative  interrogate-reply  method  such  as  is  used  in  an 
Air  Traffic  Control  Beacon  system  or  an  IFF  system.  Several  such  satellite 
based  navigation  systems  were  studied  at  Lincoln  Laboratory  (Ref.  1).  One 
further  distinguishes  such  systems  into  simplex  and  duplex  systems  depending 
on  whether  interrogate-reply  happens  sequentially  or  simultaneously  (using 
frequency  or  code  orthogonality  for  example) . 
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The  second  type  of  ranging  system  is  one  in  which  the  several  sources 
transmit  time-marked  signals  from  which  the  user  deduces  one-way  range.  If 
the  source  location  is  unknown  or  continuously  changing  as  is  the  case  for 
satellites,  the  signals  must  carry  sufficient  information  to  make  it  known. 

The  data  must  also  allow  the  time  markings  to  be  expressed  in  source-clock 
time.  The  latter  allows  the  user  to  deduce  the  time  of  transmission  of  the 
signals.  The  ranging  task  is  complicated  by  the  fact  that  the  transmitted 
time  is  expressed  in  source-clock  time  and  the  time  of  reception  in  receiver- 
clock  time.  If  however,  one  considers  the  several  source-docks  to  be  syn¬ 
chronized  (maybe  after  adding  corrections  obtained  from  the  transmitted  data) 
the  only  error  derives  from  the  user  clock.  Assuming  for  a  moment  that  all 
range-readings  are  taken  simultaneously  then  the  user  clock  error  constitutes 
just  one  more  unknown  to  be  determined  together  with  three  position  coordinates. 
A  set  of  four  ranges  allows  the  user  to  obtain  a  solution.  Non-slmultaneous 
readings  will  Introduce  additional  errors  due  to  user  clock  imperfections.  The 
GPS  C/A  and  P  modes  are  ranging  systems  of  this  second  type  and  this  is  the  only 
type  considered  further. 

2.3  Satellite  Configurations,  Frequency  Choices,  Channel  Sharing 

The  choice  of  sources  leads  to  such  diverse  ranging  systems  as  Omega, 

Loran,  Transit  and  now  NAVSTAR  GPS.  It  is  beyond  the  scope  of  this  study  to 
consider  new  satellite  configurations  or  even  the  addition  of,  for  example,  a 
geo-stationary  satellite  to  the  presently  projected  GPS  configuration.  The 
same  holds  for  frequency  choice.  No  new  frequency  choices  were  considered. 
However,  since  two  frequencies  (L^  »  1575  MHz,  I^  *  1227  MHz)  have  been  allo¬ 
cated,  either  one  or  both  could  be  considered  for  possible  waveform  change. 

Since  ranging  to  several  sources  is  necessary,  the  problem  of  channel  shar¬ 
ing  arises.  The  C/A  signals  of  GPS  are  shared  via  CDMA  with  a  separation  as 
expressed  by  the  peak  cross-correlation  of  the  PN  codes  of  21.6  dB  for  maximum 
doppler.  Within  the  class  of  CDMA  waveforms  variations  on  codelength  or  code 
generation  procedures  have  been  explored  (Ref.  2,3)  in  the  context  of  a  less 
costly  GPS  receiver  (called  Spartan) .  PN  encoding  in  a  CDMA  setup  fulfills 
simultaneously  the  role  of  providing  processing  gain  for  the  matched  filter, 
and  signal  separation. 


FDMA  is  generally  undesirable  if  each  signal  already  requires  substantial 
bandwidth  for  other  reasons.  If  one  user  receiver  must  receive  all  signals, 

FOMA  does  not  seem  amenable  to  a  simple  receiver  and  is  not  considered  further. 

TDMA,  provided  the  signal  (re-)acquisition  problem  is  kept  simple,  is 
the  solution  resulting  in  the  optimal  use  of  bandwidth  and  power.  In  the  TDMA 
systems  proposed  in  this  study  each  satellite  in  turn  is  assigned  a  fixed 
timeslot.  Guardtime  of  about  18  msec  to  avoid  signal  overlap  is  provided 
corresponding  to  the  maximum  possible  range  difference  (5600  km)  of  two  satellites 
to  a  user.  TDMA  is  ideally  suited  to  a  one-channel  sequential  receiver.  Section 
VII  in  this  report  highlights  its  benefits.  Obviously  during  each  time-slot  a 
wide  variety  of  waveforms  may  be  transmitted. 

We  have  considered  one  additional  mode  of  channel  sharing,  called  PRFDMA, 
where  each  of  the  several  signals  consists  of  an  infinite  train  of  very  short 
pulses  with  different  PRF's  (Section  6). 

2.4  TDMA  Waveform  Choices 

A  choice  of  waveform  is  based  on  a  set  of  criteria,  in  this  case  ranging 
and  communication  performance,  not  only  in  terms  of  accuracy  and  bit  error 
rates  but  also  in  terms  of  easy  signal  (re-)acquisition  and  synchronization. 

And  finally  there  is  simplicity  of  implementation  of  the  receiver. 

The  requirements  for  data  communication  are  mild;  only  50  bps  must  be 
reliably  transmitted,  and  even  so  the  same  1500  bit  message  is  repeated  about 
120  times.  Such  redundancy  could  be  exploited  to  an  extent  to  coast  through 
signal  fades. 

The  requirement  for  ranging  accuracy  is  the  most  demanding.  Ranging  is  a 
function  performed  by  radars.  All  waveforms  of  interest  have  been  analyzed  and 
synthesized  in  connection  with  radar  applications.  Waveforms  are  classified 
according  to  the  nature  of  their  ambiguity  function  (Ref.  4,5).  The  ambiguity 
function  is  the  squared  envelope  of  the  matched  filter  response  to  the  received 
signal,  often  represented  as  a  surface  with  Independent  variables  the  mismatch 
in  delay  and  doppler  between  the  signal  and  a  local  replica.  Desirable  features 
of  the  surface  are  a  sharply  defined  central  peak  for  zero  mismatch  in  delay 
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and  doppler,  if  range  and  rate  is  to  be  measured,  or  a  sharp  ridge  extending 
in  the  doppler  dimension,  if  only  range  is  to  be  measured,  and  very  low  even 
sidelobes.  At  the  same  time,  the  behaviour  of  the  surface  in  the  immediate 
vicinity  of  the  peak  or  ridge  should  have,  or  allow  easy  generation  of,  an 
indication  of  the  direction  in  which  the  peak  lies ,  for  the  purpose  of  correc¬ 
tive  action.  We  discuss  waveforms  with  these  desirable  ambiguity  function 
properties  in  mind. 


2.4.1  Constant  Carrier  Single  Pulse 


The  most  basic  ranging  waveform  is  the  constant  carrier  single  pulse 
with  basic  characteristic  a  time-bandwidth  product  of  unity.  The  ambiguity 
function  is  concentrated  in  a  unity-size  cell  in  delay-doppler  space  with  a 
negligible  portion  of  the  pulse  energy  in  the  sidelobes.  Sidelobes  are  deter¬ 
mined  by  the  exact  pulse  shape.  Ranging  accuracy  relates  to  how  well  the  peak 
of  the  ambiguity  function  can  be  located.  For  an  ideal  square  pulse  the 
ambiguity  function  in  the  delay  dimension  is  triangular*  and  would  in  principle 
result  in  perfect  accuracy.  For  a  pulse  of  practical  interest,  the  accuracy 
is  often  expressed  by  a  lower  bound  on  the  standard  deviation  of  the  ranging 
error 


ffT* 
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where  0  expresses  RMS  signal  bandwidth.  In  this  expression  0  is  the  second  cen 

2 

tral  moment  of  |S(f)|  where  S(*)  is  the  Fourier  transform  of  the  modulating 
signal;  Nq  is  one-sided  spectral  density  of  white  noise;  E  is  the  signal  energy 
The  value  of  a 1  approaches  more  closely  the  asymptotic  lower  bound  for  higher 
E/N  . 


The  above  formula  is  useful  for  sensitivity  analysis.  It  shows  that  in¬ 
creasing  bandwidth  and  E/Nq  increases  accuracy.  One  way  of  increasing  E  would 
be  to  chose  a  longer  pulse.  However,  it  can  be  shown  chat  longer  pulses  result 
in  decreased  bandwidth  0,  the  two  effects  canceling  each  other  at  least  in  the 
delay  dimension.  (Accuracy  in  the  doppler  direction,  which  is  not  considered 
here,  improves  with  longer  pulses). 


*See,  for  example,  p.  486,  488  in  Skolnik,  "Introduction  to  Radar  Systems". 


11 


If  it  is  true,  within  limits,  that  a  short  pulse  and  a  long  pulse  can  give 
equal  ranging  performance,  one  may  look  at  a  regular  train  of  short  pulses  which 
contains  the  same  energy  as  the  single  big  pulse  and  which  are  coherently 
integrated.  The  result  is  a  "bed-of-spikes"  ambiguity  function  with  improved 
delay  and  doppler  resolution  capability,  but  with  the  possibility  of  incorrectly 
relating  to  the  wrong  spike.  The  central  spike  (which  has  dimension  —  as  do 
all  spikes)  reflects  the  good  properties  of  a  single  narrow  spike  in  the  delay 
dimension  and  those  of  a  fictitious  long  pulse  (as  long  as  the  train)  in  the 
doppler  dimension.  Such  a  train  is  a  first  example,  of  waveforms  where  some  form 
of  modulation  (here  amplitude  modulation)  results  in  a  large  time-bandwidth 
product. 

Other  classes  of  so-called  pulse-compression  signals  with  large  TB-products 
have  been  developed.  They  comprise  the  waveform  classes  with  thumbtack,  sheared 
ridge  and  bed-of-spikes  ambiguity  functions. 

2.4.2  Waveforms  with  Thumbtack  Ambiguity  Function* 

In  this  class  the  large  time-bandwidth  product  TB,  resulting  in  a  resolu¬ 
tion  cell  size  of  1/TB,  is  generated  by  some  non-linear  phase  function  which 
can  be  removed  by  a  matched  filter  receiver.  Generally  Doppler  mismatch 
disturbs  this  match  and  therefore  the  effect  of  signal  compression.  In  the 
class  of  waveforms  generating  thumbtack  ambiguity  function  there  is  generally 
a  large  potential  for  interference  from  high  Doppler-delay  sidelobes. 

Consider  first  the  single-pulse  waveform  and  focus  attention  on  the  phase 
modulation.  Whether  continuously  varying  phase  variation  is  considered  (result¬ 
ing,  for  example,  in  the  subclass  of  frequency-shift  codes),  or  constant-phase 
steps  (subclasses  are  phase-shift  codes  with  further  division  into  poly  phase- 
shift  codes,  or  its  simplest  version,  the  phase-reversal  codes)  the  differences 
lie  more  in  the  realm  of  implementation  than  in  the  potential  for  ranging 
performance  as  their  ambiguity  functions  are  all  of  the  same  type,  as  long  as 
there  are  no  regularities  in,  for  example,  for  the  phase-reversal  class,  the 
sequence  of  phase- reversals.  It  is  easy  to  see,  for  example  for  the  phase- 

reversal  codes,  that  a  high  time-bandwidth  product,  TB,  means  a  long  code  (T) 
and  a  high  chipping  rate  (B) . 

*For  illustration,  see,  for  example,  Reference  4,  p.  133. 


As  a  second  case  consider  a  single  pulse  but  amplitude  modulated.  The 
amplitude  may  vary  continuously,  or  in  steps.  As  long  as  there  is  no  regularity 
in  the  way  the  amplitude  varies  a  thumbtack  ambiguity  function  results.  Even 
if  the  extreme  amplitude  variation  of  on-off  is  used,  but  the  resulting  pulses 
are  irregularly  spaced,  the  irregular  time  spacing  (sometimes  called  staggering) 
will  avoid  doppler  ambiguities  becoming  pronounced  as  they  were  for  a  regular 
pulse  train.  Instead  of  irregular  spacing,  the  pulses  could  be  pseudorandomly 
coded  in  phase  or  frequency  with  a  similar  result. 

Generally  the  thumbtack  function  is  the  "ideal"  ambiguity  function 
for  ranging  even  if  the  sidelobe  levels  are  high,  although  less  pronounced 
than  is  the  case  in,  for  example,  bed-of-spikes  type  ambiguity  functions. 

2.4.3  Waveforms  with  Sheared  Ridge  Ambiguity  Function* 

This  is  the  class  principally  of  linear  FM  waveforms  or  stepped-frequency 
approximations  to  linear  FM,  and  coherent  pulse  trains  with  linear  frequency 
shifting  from  pulse  to  pulse.  The  basic  drawback  of  this  class  of  waveforms 
is  that  no  independent  measurement  capability  exists  in  range  and  doppler. 

While  a  frequency  sweep  is  easy  to  implement,  one  needs  an  aid  to  resolve  the 
range-doppler  ambiguity.  One  such  solution  is  to  have  an  independent  doppler 
reading  (for  example,  from  a  pulse  with  no  FM  modulation)  while  another  solution 
is  to  have  two  FM  signals  each  with  different  slope.  In  our  study  this  class 
was  not  further  considered  since  such  "fixes"  would  result  in  more  complicated 
designs  than  in  the  case  with  the  other  classes  discussed  here. 

j  2.4.4  Class  of  Harmonic  Signals 

This  is  the  general  class  including  harmonic  tones  and  specially  con¬ 
structed  square  waveforms  such  as  BINOR  (Ref.  6).  The  latter  replaces  the 
sinusoids  of  the  tones  by  geometrically  related  square  waves  as  far  as  the 
generation  of  the  baseband  waveform  is  concerned.  They  may  differ  in  the 
choice  of  RF  modulation  where  BINOR  could  use  DSB  with  suppressed  or  reduced 
carrier,  while  tones  use  phase  or  frequency  modulation.  The  square  waves 
approach  is  known  not  to  be  superior  to  PN  encoding  in  performance  or  ease 
of  implementation  and  was  therefore  not  studied  further. 

*For  illustration,  see,  for  example,  Reference  4,  p.  237. 
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The  cone  ranging  approach  is  interesting  because  of  the  properties  of 
the  ambiguity  function  in  the  delay  dimension.  Accuracy  is  determined  by  the 
choice  of  the  highest  frequency.  Ambiguity  (repetition  of  the  highest  peak 
where  tones  are  cophasal)  is  determined  by  the  lowest  harmonic  tone.  In 
between  peaks  the  ambiguity  function,  around  the  highest  peak,  varies  mono- 
tonically.  This  slope  is  an  indication  of  the  direction  in  which  the  delay  has 
to  be  changed  in  order  to  achieve  the  peak.  This  is  in  strong  contrast  with 
PN-encoding  where,  outside  of  the  region  of  alignment  of  about  two  chips,  no 
indication  exists  as  to  whether  to  advance  or  retard  the  local  code  replica 
for  better  code  coincidence.  Therefore,  for  long  codelengths  and  low  signal 
levels,  as  is  the  case  for  the  present  GPS  signal,  where  matched  filtering 
is  impractical  at  present,  one  has  potentially  a  serious  code  acquisition  and 
synchronization  problem.  This  is  even  more  serious  for  the  sequential  receiver 
where  the  several  codes  are  continuously  to  be  re-synchronized. 

The  potential  of  tone  ranging  to  allow  fast  synchronization  was  promising 
enough  to  warrant  a  complete  waveform  design  effort  (Section  5). 
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3.  USE  OF  PRESENT  C/A  WAVEFORM  WITH  MORE  POWER 


In  this  section  we  summarize  the  results  of  a  study  on  the  effects  on 
receiver  complexity  if  the  present  C/A  waveform  would  be  retained  but  the 
available  power  level  would  be  increased  by  a  substantial  amount  (in  the 
range  of  3  to  15  dB) . 

3.1  Context  for  the  Study  on  Increased  Power 

A  generic  GPS  receiver  is  considered  and  its  several  functions  analyzed. 
In  the  analysis  the  following  requirements  have  been  assumed: 

The  presently  proposed  C/A  GPS  waveform  must  be  used. 

The  receiver  must  be  sequential. 

The  receiver  must  have  simultaneous  ranging  and  data 
collecting  capability. 

The  output  must  be  3D  position  and  velocity  with 
accuracies  comparable  to  those  presently  advertised 
for  GPS . 

The  emphasis  on  the  use  of  more  power  is  not  necessarily 
to  provide  more  accurate  ranging  but  simplified  circuitry, 
but  neither  will  a  highly  degraded  ranging  capability 
be  acceptable. 

The  envisioned  user  class  requires  that  the  receiver  be 
capable  of  coping  with  moderate  dynamics  (lg).  Higher 
dynamics  are  allowed  to  reduce  performance  but  not  to 
make  the  receiver  totally  fail. 

The  noise  consists  of  ambient  noise  plus  unintentional 
radar  interference  in  a  1  MHz  band  at  the  frequency. 

The  GPS  waveform  dictates  certain  receiver  functions.  This  study 
addresses  the  performance  of  these  functions  when  the  power  level  is 
increased . 
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Specifically: 


The  code  synchronization  function  is  analyzed.  Of 
concern  is  the  time  it  takes  to  synchronize.  The 
sync  function  provides  the  ranging  capability  and 
allows  carrier  despreading: 

The  carrier  acquisition  and  tracking  function  is  analyzed. 
Again  acquisition  time  is  of  concern  as  is  the  sensitivity 
of  the  phase  tracking  loop  to  user  dynamics.  The  principal 
outputs  are  doppler  reading  and  data  demodulation 
capability. 

The  data  demodulation  function  is  analyzed  in  terms  of 
error  rates  as  a  function  of  signal-to-noise  ratio. 

The  sequencing  function  is  analyzed.  The  "dwell-time" 
per  satellite  determines  user  track  update  rate  and  to 
a  certain  extent,  complexity  of  tracking  algorithm. 

Other  related  Issues  are  the  relationship  of  bandwidth 
of  the  Phase  Locked  Loop  (PLL)  to  dwell  time,  and  the 
effect  of  dwell  time  on  allowable  oscillator  quality. 

3.2  Summary  Statements  on  Impact  of  More  Power  on  GPS  Receiver 

In  this  section  we  summarize  the  conclusions  of  our  study.  The  details 
of  the  analysis  are  in  Appendix  A. 

3-2.1  Code  Synchronization 

An  increase  of  power  by  an  amount  of  10  dB  would  reduce  the  synchronization 
detection  time  on  a  given  code  phase  to  as  low  as  one  msec.  That  time  is 
particularly  suited  to  the  C/A  code  since  it  corresponds  to  a  full  code  length 
of  1023  chips.  The  corresponding  doppler  which  can  be  tolerated  is  ±  500  Hz. 

The  above  calculation  is  valid  for  a  detection  threshold  of  at  least  10  dB 
over  the  1  kHz  bandwidth.  In  normal  operation  of  a  sequential  receiver  with 
dwell  time  on  each  satellite  of  the  order  of  one  second  and  a  relatively  cheap 
Voltage  Controlled  Oscillator  (VCO)  of  10  ^  Hz/Hz,  it  is  expected  that  the 
received  carrier  frequency  at  the  end  of  a  three-second  period  will  still  be 
within  the  bandwidth,  for  a  user  with  modest  (lg)  dynamics.  It  is  also  expected 
that  the  code  phase  can  be  predicted  to  within  1  or  2  chips.  Re-acquisition 
of  the  code  then  will  only  require  at  most  a  few  code  phase  tests  and  at  worst 
one  or  two  frequency  steps.  Total  re-acquisition  time  will  then  require  only 
a  few  msec. 
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The  situation  is  entirely  different  at  cold  start  when  no  a  priori  knowl¬ 
edge  is  available  on  doppler  frequency  or  code  phase.  For  the  maximum  doppler 
uncertainty  of  4342  Hz  from  satellite  motion  and  1  kHz  from  user. motion  the 
maximum  synchronization  time  could  be  reduced  to  20  seconds  for  a  10  dB  power 
increase  above  the  present  GPS  level. 

3.2.2  Possible  Use  of  Shorter  Codes 

If  there  were  any  pressing  reason  to  hold  down  code-acquisition  time,  or 
if  still  more  power  than  an  increase  of  10  dB  were  available,  one  might  consider 
the  use  of  a  shorter  code  and/or  incoherent  integration  to  great  advantage.  But 
the  use  of  shorter  codes  will  narrow  the  margin  between  code  cross  correlation 
sidelobes  and  the  code  auto  correlation  function  peak.  That  margin  determines 
multiple  code  access  capability.  The  -23*9  dB  margin  (no  doppler  present)  for 
the  present  C/A  code  shrinks  to  -17.5  dB  for  a  shorter  127  bit  Gold  code.  This 
margin  could  become  even  less  due  to  different  doppler  affecting  the  different 
codes . 


3.2.3  Use  of  riatched  Filter  for  Code  Synchronization 

The  use  of  shorter  codes  could  bring  new  technology  to  bear  on  the  code 
synchronization  problem  such  as  the  use  of  matched  filters  in  the  form  of 
charge  coupled  devices  (CCD) .  The  CCD  matched  filter  obviates  the  need  for 
code-phase  search  and  the  delay- locked  loop,  and  realizes  the  full  correlation 
processing  gain.  At  present,  however,  the  applications  of  the  programmable 
CCD  correlator  in  advanced  engineering  developments,  are  still  limited  to 
codelengths  of  up  to  256  chips. 

3.2.4  Increased  Ranging  Accuracy 

The  code  tracking  accuracy  in  terms  of  the  chip  splitting  capability 
increases  roughly  as  the  square  root  of  the  signal  to  noise.  For  example,  a 
10  dB  Increase  in  power  roughly  results  in  an  increased  chip  splitting  capa¬ 
bility  of  a  factor  of  3. 

3.2.5  Increased  Performance  of  Carrier  PLL 

A  trade-off  study  was  performed  on  increase  in  performance  of  the  carrier 
PLL  as  a  function  of  increased  power.  The  presently  considered  PLL  works  with  an 
adjustable  bandwidth  to  allow  rapid  acquisition  on  the  one  hand  and  sufficiently 
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accurate  phase  tracking  on  the  other  hand.  The  degree  to  which  bandwidth  can 
be  narrowed  depends  to  a  large  extent  on  the  dynamics  of  the  user.  Sudden 
accelerations  translate  into  accelerations  of  phase.  The  transient  phase 
behavior  and  the  steady-state  loop  errors  then  depend  on  the  order  of  the  PLL 
and  in  particular  on  the  filter  in  the  PLL.  The  result  of  increased  SNR  is 
that  the  phase  error  due  to  dynamical  effects  dominates  the  noise-induced  error 
more  and  more.  The  optimal  loop  bandwidth  becomes  larger  with  increased  SNR. 
Figure  3.1  illustrates,  for  signal-to-noise  ratios  at  the  input  of  the  PLL  of  25 
to  45  dB,  what  the  optimal  choice  of  bandwidth  would  be  when  the  phase  error 
components  are  weighted  (weight  2  for  steady-state  error  resulting  from  acceler¬ 
ation  versus  weight  1  for  standard  deviation  of  the  noise  induced  phase  error) . 
The  analysis  applies  to  the  second  order  PLL.  Larger  optimal  bandwidths  here 
mean  for  example  that  the  receiver  becomes  less  sensitive  to  carrier  frequency 
offsets  (either  from  motion  or  from  clock  errors)  or,  for  a  given  offset,  can 
phase- lock  in  a  shorter  time.  Thi3  is  important  for  the  sequential  receiver 
where  one  switches  to  a  different  signal  at  a  rate  of  maybe  once,  and  possibly 
many  times,  a  second.  Initial  acquisition  will  require,  in  a  cold-start  mode, 
fewer  frequency  bands  to  be  searched. 

3.2.6  Bit  Error  Rates 

An  increase  in  signal-to-noise  ratio  contributes  directly  and  indirectly 
to  lowering  the  bit  error  rate.  First  it  increases  the  detection  energy 
available  after  a  given  integration  time.  And  secondly,  the  use  of  a  wider 
bandwidth  in  the  PLL  gives  substantial  protection  against  cycle  slip  in  a 
dynamic  user  environment,  thus  lowering  the  occurrence  of  bit  errors. 

3.2.7  Shorter  Dwel 1-Time  per  Satellite 

In  a  sequential  receiver  higher  SNR  means  a  quicker  code  and  carrier 
acquisition.  This  feature  allows  shorter  dwell  time  per  satellite  to  be  more 
efficient  in  terms  of  collecting  ranging  information  (but  not  necessarily  for 
data  collection) .  Shorter  dwell  time  will  result  in  a  higher  ranging  and 
position  update  rate.  For  a  given  requirement  on  the  accuracy  of  the  smoothed 
user  position  a  higher  update  rate  allows  individual  range  measurements  to  be 
less  accurate.  This  may  allow  a  simpler  circuit  for  delay  lock  loop  tracking 
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PIG.  3-1  MINIMUM  PHASE  ERROR  AND  OPTIMAL  LOOP  BANDWIDTH  AS  A  FUNCTION  OF  C 


to  be  sufficient.  In  the  military  GPS  receiver  implementations,  the  delay  lock 
discrimination  curve  is  usually  obtained  from  the  output  of  two  or  more  corre¬ 
lator  channels.  One  more  correlator  channel  is  used  to  perform  the  function  of 
stripping  the  code  from  the  signal  and  to  provide  the  input  to  the  Costas  loop 
for  carrier  cracking  and  data  demodulation.  Simpler  Delay  Locked  Loop  (DLL) 
error  generating  methods,  such  as  the  use  of  tau-dither  methods,  may  eliminate 
the  use  of  several  correlator  channels. 

3.2.8  Possibility  of  Using  Less  Costly  VCO 

In  continuation  of  the  argument  made  under  point  3.2.7,  it  can  be  stated 
that  if  shorter  dwell  times  can  be  implemented,  the  time  between  "visits"  to 
the  same  satellite  are  shorter.  Less  frequency  drift,  will  then  accumulate 
with  respect  to  the  carrier  frequency  of  a  particular  satellite  while  other 
satellite  signals  are  collected.  This  allows  the  use  of  a  less  stable  and 
consequently  cheaper  oscillator.  For  the  GA  application  our  designs  have 
aimed  at  1  g  user  dynamics,  a  10  ^  Hr/Hz  oscillator  and  position  update  rates 
from  1  to  4  seconds  (dwell  times  .25  to  1  sec.). 

These  eight  points  summarize  our  findings  on  the  effect  of  a  GPS  C/A 
waveform  with  more  power. 
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4.  A  TDMA  PULSED  GPS  WAVEFORM 


The  advantages  of  TDMA  waveforms  as  a  class  are  rather  compelling  for  a 
low-cost  GPS  receiver,  and  they  are  discussed  in  some  detail  in  Section  7. 

In  this  section  we  Introduce  a  particular  TDMA  example,  namely  a  pulsed  system. 
The  fundamental  tradeoff  between  the  extremes  of  pulsed  and  CW  waveforms  lies 
between  the  cost  and  difficulty  of  generating  high  power  pulses  on  board  the 
satellites  and  the  advantages  of  simplicity  and  performance  gained  for  the 
receiver.  In  this  report  we  emphasize  the  advantages  for  the  receiver,  which, 
in  our  opinion,  are  sufficiently  great  to  balance  the  possible  difficulties 
associated  with  high  transmitter  power.  A  brief  survey  has  indicated  to  us 
that  suitable,  space-qualified  power  generation  hardware  is  or  soon  will  be 
available,  and  the  interoperability  of  such  a  pulsed  transmitter  with  the 
other  functions  of  the  GPS  satellite  hardware  does  not  seem  to  present  insur¬ 
mountable  difficulties,  but  this  last  point  has  not  been  analyzed  in  detail 
in  this  study. 

Our  main  intent  was  to  complete  the  design  of  a  pulsed  GPS  waveform  and 
derive  performance  from  it  which  is  considered  representative  for  the  class  of 
pulsed  systems.  Although  an  attempt  was  made  to  select  optimum  choices  many 
refinements  are  possible. 

4.1  Waveform  Description 

In  the  system  described  here,  the  satellites  transmit  part  of  their  message 
using  short  bursts  of  high  energy  pulses .  The  message  portion  is  preceded  by  a 
set  of  identical  pulses  from  which  ranging  information  is  derived.  The  L-band 
carrier  is  spread  by  a  PN  code  (BPSK  modulation)  at  a  4  MHz  chipping  rate.  The 
data. is  PP  (Pulse  Position)  modulated  on  the  spread  carrier.  The  message  con¬ 
tent  is  chosen  to  be  the  same  as  for  the  presently  proposed  GPS  system. 

4.1.1  Power  Level 

The  first  consideration  is  the  needed  power  level.  The  signal  carrier 
power  to  noise  power  density  ratio  (C/Nq)  at  the  receiver  front  end  is  practi¬ 
cally  determined  by  the  required  detection  energy  (corresponding  to  chosen 
values  for  probabilities  of  detection  and  false  alarm)  and  the  frequency  offset 
one  must  cope  with  (which  fixes  integration  time). 


For  a  chosen  value  of  .97  probability  of  detection  and  10  ^  probability 
of  false  alarm,  based  on  coherent  integration,  the  required  signal  energy  to 
noise  power  density  ratio  Eq/No  is  14  dB.  To  simplify  the  particular  TDMA 
pulsed  design  under  consideration,  we  wanted  to  operate  at  a  fixed  receiver 
frequency,  independently  of  the  doppler  present,  or  the  clock  error.  To  main¬ 
tain  an  Eq/n0  of  14  dB  under  the  worst  frequency  offset  of  5.5  kHz  with  the 
above  probabilities  the  Integration  time  was  chosen  to  be  64  ysec  and  the  power 
ratio  C/Nq  to  be  58  dBHz.  The  Eq/No  for  no  frequency  error  is  then  16  dB.  For 
a  receiver  noise  figure  of  5  dB  the  required  received  power  is  -141  dBw. 

Choosing  a  receiver  antenna  gain  of  0  dB  and  a  satellite  antenna  gain  of 
11.5  dB,  a  path  loss  of  184.5  dB  at  the  longest  user-satellite  distance  at 
5°  elevation,  and  an  atmospheric  loss  of  1  dB,  the  required  transmitted  peak 
power  would  be  2  kw.  Table  4-1  summarizes  the  power  budget. 

4.1.2  Signal  Format 

We  chose  a  strict  TDMA  format,  allowing  only  the  two  satellites  in  anti¬ 
podal  positions  to  transmit  simultaneously.  The  total  timeframe,  of  as  yet 
unspecified  length,  will  therefore  be  divided  into  12  slots.  Since  a  time 
measurement  will  be  obtained  from  at  least  four  satellites  per  frame,  trans¬ 
mitting  in  designated  slots,  the  position  update  rate  is  determined  by  the 
frame  length.  Update  rates  should  be  higher  than  once  per  several  seconds  for 
adequate  tracking  performance.  An  argument  can  be  made  that  higher  update 
rates  allow  good  tracking  through  maneuvers  or  allow  the  tracking  filters  to 
be  substantially  simpler.  Most  importantly  a  high  update  rate  or  short  frame 
time  avoids  the  possibility  of  the  short-term  clock  drift  accumulating  any 
substantial  timing  error.  The  issue  of  clock  stability  in  relation  to  these 
receivers  is  the  subject  of  Section  9.  We  therefore  opt  for  short  frames. 

In  the  design  we  chose  a  one  half  second  frame  length.  In  a  strict  TDMA 
system  the  timeslot  per  satellite  pair  is  then  41.6  msec.  To  avoid  overlap 
between  the  transmissions  of  different  satellites  (closest  at  zenith  and  far¬ 
thest  at  lowest  elevation)  a  "minimum"  guard  time  of  17  msec  is  needed.  We 
expanded  this  to  19.134  msec  after  other  time  parameters  were  chosen.  The  useful 
transmission  time  is  divided  into  preamble  and  data  transmission  as  shown  in 
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Transmitted  Peak  Power 

+  33  dBw  (2  kw) 

Satellite  Antenna  Gain 

+  11.5  dB 

Path  Loss  on  frequency  L1  on 

Longest  Path  (25824  km) 

-  184.5  dB 

Atmospheric  Absorption 

-1  dB 

Receiver  Antenna  Gain 

0  dB 

Received  Peak  Power 

141  below  lw  or  -141  dBw 

Receiver  Noise  Figure  of  5  dB 

199  dB  below  lw  per  Hz 
or  -199  dB- joule 

Ratio  of  Received  Peak  Power  to 
Noise  Power  Density 

58  dBHz 

Per  Pulse  of  64  ysec 

-42  dB  sec 

Ratio  of  Received  Signal  Energy  to 
Noise  Power  Density 

14  to  16  dB  (up  to  5.5  kHz 
frequency  offset  allowed). 

■ST 


The  preamble  consists  of  eight  pulses,  allowing  synchronization  with  the 
PN  code  and  accurate  timing  (ranging) .  The  message  portion  consists  of  twenty 
pulses,  each  pulse  conveying  four  data  bits  (see  Section  4.1.5  on  modulation). 

4.1.3  Waveform 

All  pulses  are  64  psec  long  and  identical.  The  1^  carrier  is  BPSK  modu-  . 
lated  by  a  code  256  chips  long,  at  a  4  MHz  chipping  rate.  The  higher  chipping 
rate  was  chosen  over  the  1  MHz  rate  for  the  present  C/A  GPS  signal  in  order  to 
reduce  the  receiver's  sensitivity  to  multipath  by  a  factor  of  four.  For  a  4  MHz 
rate  the  altitude  above  which  no  multipath  effects  from  a  flat  surface  are  felt 
reduces  to  about  2,000  ft.  for  a  minimum  elevation  angle  of  the  satellite  of  5°. 

Since  the  signals  of  the  12  satellite  pairs  are  strictly  time  multiplexed, 
a  common  single  code  can  be  used. 

4.1.4  Ranging  and  Synchronization 

As  shown  before,  the  preamble  consists  of  8  pulses.  We  consider  here  the 
issue  of  detection,  synchronization  and  ranging. 

Table  4-1  showed  the  SNR  per  pulse  at  the  output  of  the  matched  filter  to 
be  14  dB  with  2  dB  margin.  This  establishes  a  probability  of  detection  (inco¬ 
herent  phase)  of  PD  ■  97%  and  a  probability  of  false  alarm  of  10 

The  2  dB  margin  is  used  against  the  drop  in  SNR  due  to  doppler  mismatch. 

It  can  be  easily  shown  that  for  a  worst  doppler  of  ±  5500  Hz  the  power  falloff 
of  the  autocorrelation  function  is  1.8  dB  for  a  64  psec  integration  time. 

The  range  accuracy  in  terms  of  chip-splitting  capability  is 

a  -  A  —  1  -  -  0.1A 

2/E7iT 

o 

where  is  the  standard  deviation  of  the  TOA  error  from  receiver  noise  only. 
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The  sampling  noise,  at  8  MHz  sampling  rate,  adds 

a  =  A  — - — 

s  2/12 

The  basic  accuracy  per  pulse,  for  A  ■  250  nsec  is  then 


/  2  2“ 

a  +  a  -  44  nsec. 

T  S 

However,  ranging  will  be  based  on  the  timing  of  8  pulses  resulting  in 


a 


R 


1 

/r 


15.5  nsec  or  ft  » 


which  is  the  basic  ranging  accuracy.  Obviously  other  error  sources  add  to 


this  error.  One  is  the  errors  n  satellite  position  which  is  implicitly  trans¬ 
mitted  in  the  data  and  which  can  be  calculated  by  the  user .  These  "  ephemeris 
errors"  are  estimated  to  be  smaller  than  5  feet.  The  main  reason  for  keeping 
these  errors  so  small  is  that  the  identical  data  are  used  in  the  P-mode  with 
basic  ranging  accuracies  claimed  better  than  5  ft  and  all  other  errors  must 
therefore  be  kept  commensurately  small.  Errors  from  tropospheric  delays  and 
satellite  clocks  are  of  the  order  of  a  few  feet.  Errors  from  multipath  can  be 
a  substantial  fraction  of  a  chip.  For  this  reason  the  higher  chipping  rate  of 
4  MHz  was  chosen.  For  this  reason  the  altitude  floor  above  which  a  user  would 
not  be  affected  by  multipath  is  reduced  from  8000  ft  to  about  2200  ft.  The 
principal  potential  error  contribution  is  from  uncorrected  ionospheric  delay. 
Since  a  simple  sequential  GPS  receiver  probably  does  not  use  the  dual  frequency 
correction  method  one  must  rely  on  a  model.  Several  such  models  have  been 
proposed.  It  was  estimated  that  with  the  use  of  a  model  statistically  50%  of 
the  error  remains.  One  redeeming  feature  is  that  the  error  is  strongly  corre¬ 
lated  in  time,  thus  acting  more  like  a  temporary  bias  rather  than  a  random  error. 
The  resulting  user  position  will  be  slightly  in  error  but  in  a  consistent  way. 
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The  model  likely  to  be  used  is  similar  to  the  General  Dynamics  "Semi-Static 
Crossline  Ratio  Technique."  The  technique  is  based  on  two  tables  stored  in  the 
user’s  processor  and  which  correspond  to  the  two  seasonal  extremes.  Each  table 
gives  vertical  group  delay  for  5  degree  increments  around  the  magnetic  equator- 
For  every  5  degrees  increment  in  latitude  away  from  the  equator  it  provides  the 
ratio  of  the  vertical  delay  as  related  to  the  corresponding  delay  at  the  mag¬ 
netic  equator.  Satellite  data  provide  three  coefficients  for  seasonal  inter¬ 
polation.  The  delay  at  any  angle  different  from  vertical  can  then  be  calculated. 

4-1.5  Data  Modulation 

With  a  waveform  based  upon  a  single  coded  pulse,  it  is  natural  to  think 
of  pulse-position  modulation  (PPM)  for  data  transmission  because  of  its  economy 
of  power.  In  other  words,  a  given  bit  error  rate  (BER)  can  be  achieved  with 
less  signal-to-noise  ratio  than  higher  duty  cycle  waveforms.  The  message  block 
is  at  most  22.5  msec  in  length,  since  the  preamble  requires  2.048  msec  and  the 
minimum  dead  time  is  17.09  msec.  This  message  block  is  roughly  350  pulse 
lengths  in  duration,  and  the  only  problem  is  to  divide  this  interval  into 

£ 

symbols,  each  of  which  conveys  K  bits  by  placing  each  pulse  in  one  of  M  *  2 
positions.  Then  MN  -  350,  where  N  is  the  number  of  symbols  per  message  block. 

Pulse  position  modulation  is  a  special  case  of  M-ary  signaling  with  ortho¬ 
gonal  waveforms.  In  the  incoherent  case  (appropriate  here)  the  probability  of 
a  symbol  error  (PSE)  is  given  by  the  expression 

2  2 
-a  M  a 

PSE(M)  *  I  e  2  2  ,  (M)(-l)m  e2m  > 

M  m*2  m 

where  M  is  the  number  of  pulse  positions  and 
2  E 

—  ■  _JL  »  energy  per  pulse 

2  Nq  noise  power  density 

A  symbol  error  causes  a  number  of  bit  errors  and  the  probability  that  an 
arbitrary  bit  is  in  error,  is  given  by 

BER  “  i  M=f"  PSE  • 
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Allowing  a  2  dB  margin,  our  signal-to-noise  ratio  is  14  dB,  and  for  this  value 

of  E  /N  the  BER  varies  with  M  as  shown  in  the  following  table: 
p  o  ^ 

i  , 

K  M  ;  BER 

2  4  3.48  x  10"6 

3  8  6.86  x  10~6 

4  16  1.33  x  10-5  . 

The  larger  the  value  of  M,  the  lower  the  duty  cycle  of  the  waveform,  but  it  is  de¬ 
sirable  to  keep  the  BER  on  the  order  of  10  ^  or  lower,  hence  M“16  is  the  logical 
choice.  Actually,  since  bit  errors  occur  in  groups  with  PPM,  the  BER  as  defined 
above  is  somewhat  misleading.  A  better  criterion  would  be  expressed  in  terms  of 
mean  time  between  bit  errors,  i.e.,  the  expected  length  of  a  run  of  correct  bits. 

With  M*16  and  Ep/NQ  ■  14  dB,  the  symbol  error  probability  is  PSE  *  2.49  x  10 

and  it  is  easy  to  see  that  this  yields  the  same  error-free  run  length  as  an  equi¬ 
valent  code  with  independent  bit  errors  and  a  BER  of  6.23  x  10~^.  Hence  four 
bits  per  symbol  is  an  adequate  choice. 

Since  a  symbol  is  now  1.024  msec  long,  the  message  block  can  hold  no  more 
than  21  symbols,  and  a  raasonable  value  is  20  symbols,  which  makes  the  message 
block  20.48  msec  long  and  leaves  19.14  msec  for  the  guard  time.  The  design 

then  provides  80  bits  per  frame,  which  can  be  broken  into  5  bits  for  satellite 

ID  and  75  bits  for  data,  for  an  overall  data  rate  of  150  bps.  Table  4-2.  sum¬ 
marizes  the  data  demodulation  aspect  of  the  pulsed  design.  Figure  4-2  illustrates 
for  the  chosen  16-ary  PP  modulation  how  the  symbol  error  rates  depend  on  SNR. 

The  system  design  provides  14  dB  SNR  with  a  2  dB  margin.  Even  if  the  margin 
is  exhausted  the  symbol  error  performance  is  quite  satisfactory. 
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TABLE  4-2.  DATA  DEMODULATION  SUMMARY 


Data  Message  Length: 

1500  bits 

Data  Rate  : 

80  bits/frame  of  .5  sec  with 

- 

5  bits  for  ID  and 

- 

75  bits  of  data 

At  14  dB  SNR 

probability 
of  symbol  error  : 

2.49  10~5 

Average  length  of 
error-free  run: 

4 

1  4.02  10  symbols 

\  -  1.61  105  bits 

1  m  107  messages 

*  17.8  minutes 

Bit  error  rate  for  independent 
errors: 

6.23  10  ®  per  sec 

I  . 
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AVERAGE  MESSAGE  RUNLENGTH  WITHOUT  ERROR  VERSUS  SNR  FOR 


4.2  Receiver  Design 


The  following  section  summarizes  the  results  of  a  detailed  look  at  the 
implementation  aspects  such  as  choice  of  hardware,  size  and  to  a  limited 
extent  cost  for  the  TDMA  waveform  using  PN  encoded  pulses. 

The  hardware  implementation  sketched  out  is  only  one  of  many  approaches 
and  is  not  intended  as  the  model  of  future  hardware  but  rather  as  a  representa¬ 
tion  of  one  possibility  for  purposes  of  sizing  (and  to  an  extent  costing  as  well) . 

4.2.1  Receiver  Hardware  Design 

The  overall  receiver  block  diagram  is  shown  in  Fig.  4-3.  The  receiving 
system  consists  of  three  sections;  the  IF/Demodulator  section,  the  Special- 
Purpose  Logic,  and  the  Microcomputer.  The  IF/Demodulator  filters  and  ampli¬ 
fies  incoming  RF  signals,  heterodynes  the  desired  signal  to  baseband,  demodu¬ 
lates  the  PN  coding  in  a  correlator-pair  which  provides  integration  over  the 
pulse  duration,  and  converts  the  resultant  correlator  output  into  a  digital 
output  sample.  The  correlator  output  is  also  fed  to  the  sync  detection  circuitry. 

The  special-purpose  logic  provides  three  functions:  sync  detection,  time- 
of-arrival  determination,  and  data  demodulation.  This  information  is  fed  to 
the  microcomputer. 

The  microcomputer  performs  the  data  verification,  selection  of  range 
samples  to  be  used,  and  determination  of  position.  As  noted  in  the  intro¬ 
duction,  this  study  did  not  consider  the  computational  aspects  of  position 
determination;  that  work  remains  to  be  done  at  a  later  date. 

A  block  diagram  of  the  IF /Demodulation  circuitry  prior  to  the  correlators 
is  shown  in  Fig.  4-4.  The  design  approach  is  straightforward.  The  most 
difficult  item  is  probably  the  relatively  narrow-band  filter  used  on  the 
front  end  to  reject  interfering  signals. 

The  baseband  correlator,  combining,  and  analog-to-digital  conversion 

circuitry  is  essentially  the  same  as  that  used  in  the  development  of  the 
*  / 

JRSVC  radio  .  The  direct  applicability  of  that  work  would  reduce  the  effort 
required  to  configure  the  GPS  receiver. 

A  block  diagram  of  the  special-purpose  logic  is  shown  in  Fig.  4-5. 

* 

Work  presently  underway  at  Lincoln  Laboratory. 
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FIG.  4-3  OVERALL  BLOCK  DIAGRAM 


FIG.  4-5  SPECIAL  PURPOSE  LOGIC 


The  analog  output  of  the  correlator  combining  circuitry  is  fed  into  the 
sync  'delay  circuit,  wher.*  it  is  thresholded  and  shifted  into  a  shift  register. 

As  signals  propagate  through  the  shift  register,  a  point  in  time  will  occur 
when  sync  pulses  appear  at  all  the  taps  of  the  register  simultaneously,  and 
this  is  the  Instant  at  which  receipt  of  sync  is  recognized.  The  summed  out* 
put  of  the  register  taps  is  compared  to  one  of  two  possible  values,  and  an 
output  declared  when  the  threshold  is  exceeded.  The  reason  for  the  use  of 
two  threshold  values  is  to  allow  signals  to  be  acquired  only  when  all  eight 
sync  pulses  are  detected,  but  to  remain  in  track  with  less  than  eight  detections. 
This  reduces  the  probability  that  the  signal  will  fall  from  track  merely  because 
one  sync  pulse  was  missed. 

The  detected  sync  starts  the  data  demodulation  counter  which  controls 
the  demodulation  of  the  data,  and  also  loads  the  time-of -arrival  register 
which  stores  the  time  at  which  the  sync  arrived.  The  data  demodulation 
counter  causes  strobe  signals  to  be  generated  which  sample  the  A/D  output 
at  the  sixteen  possible  data  correlation  times  and  determine  the  largest, 

or  most  probable,  position  of  the  data  pulse.  This  position,  which  yields 

* 

the  data,  is  stored  and  subsequently  read  out  to  the  computer. 

During  initial  acquisition  the  first  sync  signal  received  starts  data 
demodulation  and  the  ID  signal  recovered  from  the  data  is  loaded  into  the 
segment  counter  to  initialize  the  receiver  with  respect  to  the  transmit 
pattern.  IDs  of  subsequent  signals  should  then  agree  with  the  segment 
counter  reading  at  the  time  they  are  received;  if  not,  a  misalignment  is 
present,  and  acquisition  is  reinitiated.  Refinements  of  the  alignment  veri¬ 
fication  can  readily  be  incorporated  if  desired.  . 

Data  is  fed  out  from  the  special-purpose  logic  to  the  microcomputer  at 
the  end  of  each  data  symbol  time.  The  microcomputer  performs  data  manipula¬ 
tion  and  verification,  as  well  as  processing  of  time-of-arrival  information 
to  determine  position. 


♦Detailed  diagrams  on  sync  processing  logic  resulting  in  logic  cards  labeled 
numbers  1  to  3  are  included  as  Appendix  B. 


35 


* 


4.2.2  Receiver  Packaging  and  Sizing 

'In  order  to  arrive  at  some  estimate  of  receiver  complexity  and  size  it 
was  necessary  to  generate  a  first  pass  at  the  design  and  then  estimate  what 
it  implies  in  terms  of  a  receiver  package.  The  receiver  front  end  and  special- 
purpose  logic  were  considered  in  some  detail,  and  drawings  generated  are  in¬ 
cluded  in  the  Appendix.  Other  portions  were  estimated  based  on  similar  work 
on  the  JRSVC  radio  or  other  projects  at  Lincoln  Laboratory. 

The  receiver  packaging  as  presently  conceived  is  shown  in  Fig.  4-6.  The 
"receiving"  portion  of  the  package,  containing  circuitry  prior  to  the  micro¬ 
computer  is  approximately  4.75x10x3“  142.5  cubic  inches.  In  the  overall 
package  as  shown  provision  was  made  for  three  microcomputer  boards.  Depending 
on  the  speed  of  computation  required,  and  the  amount  of  memory  required,  this 

may  not  be  sufficient.  If  not,  the  package  would  obviously  grow  beyond  the 

3 

4.75  x  10  x  4.25  size  (203  in.  )  shown  in  Fig.  4~6. 

A  rule  of  thumb  for  weight  of  equipment  of  this  type  is  25  cubic  inches 
per  pound,  which  would  lead  to  a  receiver  weight  estimate  of  approximately 
6  pounds. 

In  summary,  this  initial  broad-brush  look  at  the  design  of  a  receiver  for 
use  with  the  proposed  Pulsed  GPS  system  indicates  that  adoption  of  a  system  of 
this  type  could  lead  to  implementation  of  navigation  receivers  for  the  civilian 
market  at  moderate  prices. 
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FIG.  4-6  PACKAGING  OF  TDMA  PULSED  GPS  RECEIVER 


4.3  Conclusions  on  the  TDMA  Pulsed  Waveform  Design 

We  have  examined,  out  of  the  TDMA  waveforms,  the  class  of  pulsed  signals 
and  composed  a  signal  structure  which,  together  with  the  emerging  CCD  techno¬ 
logy,  has  led  us  to  a  design  resulting  in  a  simple  receiver.  The  main  advantages 
for  the  hardware  are  the  absence  of  any  DLL  or  PLL.  Also  the  pulses  are  short 
enough  to  be  insensitive  to  several  kHz  of  frequency  offset  either  due  to 
relative  user-satellite  motion  (Doppler)  or  oscillator  instability. 

Certain  features  such  as  the  4  MHz  signal  bandwidth,  which  is  beneficial  in 
obtaining  high  ranging  accuracy,  were  chosen  more  for  protection  against  multi- 
path.  The  ceiling  of  approximately  8600  ft.  above  which  there  is  no  marked 
effect  from  ground  related  multipath  for  the  GPS  C/A  signal  is  reduced  to 
2150  ft.  So  while  we  increased  the  bandwidth  by  a  factor  of  four,  it  is  still 
less  than  half  that  of  that  allotted  to  the  GPS  P-mode  on  the  same  frequency.  At 
the  same  time,  an  average  transmitted  power  of  about  one  third  that  of  the  C/A 
mode  is  needed  to  provide  very  satisfactory  performance  for  ranging  and  data 
transmission.  The  data  rate  provided  in  the  design  was  three  times  the  present 
C/A  GPS  signal  data  rate  mainly  to  reduce  the  time-to-first-fix  (TTFF) . 

The  PN-code  chosen  was  only  255  chips  long  compared  to  1023  for  the  C/A 
GPS  signal.  The  latter  has  interference  from  other  GPS  satellite  signals  to 
deal  with  (cross  correlation  sidelobes)  while  in  the  candidate  design  there  is 
only  one  signal  at  the  time.  The  number  of  chips  in  a  pulse,  whose  length  was 
determined  by  required  energy  per  pulse  for  detection  and  the  desire  to  be 
insensitive  to  the  maximum  occurring  doppler,  was  then  more  a  question  of 
accurately  determining  TOA  of  the  pulse  in  noise.  The  above  limitation  due  to 
doppler  and  a  low  14  dB  SNR  per  pulse  has  been  improved  upon  in  the  design  by 
using  a  pulse  burst  (here  8)  to  synchronize  for  data  demodulation  and  obtain¬ 
ing  an  Improved  measurement  of  pseudorange. 

Finally,  pseudorange  readings  on  all  visible  satellites  are  obtained 
within  .5  seconds.  This  provides  the  capability  of  varying  user  tracker 
update  rates,  depending  on  the  computational  capabilities  provided.  For 
example,  a  one  second  update  rate  may  be  desirable  during  maneuvers,  while  a 
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fout  second  update  interval  may  be  quite  acceptable  in  normal  flight.  This 
capability  in  turn  may  reduce  the  complexity  of  the  user  tracker.  The 
.5  second  pseudorange  update  rate  allows  several  such  measurements  to  be 
made  per  user  tracker  update,  providing  the  capability  of  extrapolating  (or 
interpolating)  the  measurements  of  the  several  satellites  to  a  common  time. 

This  feature  compensates  for  the  loss  of  pseudorange -rate  measurement  available 
to  the  C/A  GPS  receiver  with  a  phase- locked  loop.  The  interpolation  process 
is  discussed  in  detail  in  Section  9. 


TABLE  4-3.  SUMMARY  OF  FEATURES  AND  PERFORMANCE  OF  TDMA  PULSED  GPS  DESIGN 


Power  Levels 

Transmitted  Average  Power  7.2  w 

Transmitted  Peak  Power  2  kw 

Transmission  Format  8  pulse  burst  over  2.048  msec.  (25%  short  term  duty 

cycle),  20  pulses  spread  over  20.48  msec. (6.25% 
short  term  duty  cycle) . 

Deadtime  of  477  msec.  (.35%  overall  duty  cycle). 


Power  at  Receiving  Antenna  -141.25  dBw 
Received  C/Nq(5  dB  noise  figure)  56  dB-Hz 
E/N  per  pulse  14  dB 

O' 


Signal  Design  Features 

Bandwidth  at  baseband 

Ranging  Capability  Provided  by 

Data  Carrying  Method 

Data  Rate 

Data  Repetition  Period  for  user 

4  MHz 

sync  on  PN  encoded  pulse  burst 

16-ary  PP-Modulated  pulses 

150  bps 

10  sec 

Receiver  Features 

Ranging 

CCD  matched  filter  and  special  sync  logic 

Carrier  Tracking  Loop 

Absent 

Data  Demodulation 

PPM  demodulation  logic 

Start-up  Circuitry 

Absent 

Cold  Start  Capability 

Superior  to  capability  of  C/A  waveform 

Range  Update  Rate 

2  Hz  simultaneously  for  all  visible 

t 

satellites 

Measurement  Type 

Range 

Sensitivity  to  Multipath 

None  above  2200  ft  originating  from 
ground  reflections 

5.  A  TDMA  TONE  RANGING  GPS  WAVEFORM 

As  pointed  out  in  Section  2,  Tone  Ranging,  based  on  the  properties  of  the 
autocorrelation  function,  has  the  potential  to  be  a  very  good  ranging  waveform 
if  sufficient  range  ambiguity  resolution  is  present.  The  latter  is  provided 
by  using  several  tones.  The  tones  can  be  chosen  to  be  clustered  to  maximize 
bandwidth  utilization.  Tone  ranging  is  vulnerable  to  jamming.  However,  in  this 
study,  jamming  resistance  or  security  were  not  considered. 

Any  alternative  to  the  present  continuous  PN  encoded  waveform  must  provide 
ranging  and  data  carrying  capabilities.  In  this  section  the  ranging  capability 
is  provided  by  a  set  of  tones  frequency  modulating  a  carrier.  Data  are  effi¬ 
ciently  transmitted  in  separate  time  intervals,  BPSK  modulating  the  same 
carrier. 

Tone  ranging  is  an  established  technique  (Ref.  6).  A  well-known  system  is 
GRARR  (Goddard  Range  and  Range  Rate  System)  built  in  the  mid-sixties.  The  sys¬ 
tem  used  a  set  of  harmonic  sidetones  for  tracking  spacecraft  in  near-earth 
orbits  and  a  combination  of  tones  and  PN  coded  signals  for  cislunar  and  trans- 
lunar  distances.  Studies  on  tone  ranging  have  also  been  performed  in  support 
of  the  AEROSAT  program  (Ref.  7). 

Several  things  have  been  learned  from  these  studies.  We  incorporate  in 
our  design  some  features  without  further  motivation.  The  first  one  is  the  use 
of  clustered  tones  instead  of  the  fundamental  tones  needed  for  range  ambiguity 
resolution.  The  second  feature  is  the  method  for  carrying  data  on  a  separate 
subcarrier  or  on  the  basic  carrier  itself  in  a  time  multiplexed  mode.  A  third 
one  is  the  choice  of  modulation  (PM  or  FM) .  The  constant  signal  envelope,  in 
that  case,  allows  the  output  amplifiers  to  be  operating  at  their  saturated 
level  resulting  in  maximum  output  and  efficiency.  Although  exponential 
modulation  is  used  at  transmission,  no  complementary  demodulation  can  be 
used.  The  reason  is  that  for  the  typical  wide  bandwidth  of  such  modulations 
the  received  signal  to  noise  levels  (below  0  dB)  do  not  allow  regular  FM  or  PM 
demodulation  to  be  effective.  Instead,  a  direct  filtering  of  the  sidetones, 
either  at  IF  or  even  better  (3  dB),  at  baseband,  is  used. 
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The  chief  difference  between  the  GPS  tone  ranging  problem  and  previous 
tone  ranging  systems  is  that  not  one  but  twenty-four  radiating  sources  are 
involved.  Each  satellite  must  use  basically  the  same  set  of  fundamental  tones. 
If  all  satellites  use  the  same  carrier  frequency  as  well,  a  serious  inter¬ 
ference  problem  will  exist.  Indeed  the  FM  generated  RF  spectrum  consisting  of 
the  tones,  harmonics  and  intertone  harmonics  is  very  wide,  the  width  being 
determined  by  the  modulation  index.  The  only  way  continuous  signals  would  not 
interfere  would  be  when  the  corresponding  satellites  are  not  simultaneously 
visible  to  the  user  or  the  difference  in  doppler  is  sufficient  to  separate  the 
spectral  lines  for  the  receiver. 

A  solution  to  this  problem  is  to  use  a  TDMA  waveform  where  only  two 
satellites  (in  antipodal  positions)  transmit  simultaneously.  Each  pair  of 
satellites  then  transmits  in  one  of  twelve  assigned  timeslots.  The  twelve 
timeslots  make  up  the  basic  cycle  time  which  will  determine  the  position  up¬ 
date  rate  for  the  user. 

Another  solution  would  be  to  have  each  pair  of  satellites  transmit  on  a 
different  carrier  (FDMA) .  Since  the  useful  spectrum  consists  only  of  a  few 
lines  it  may  be  possible  to  "pack"  the  spectrum  with  12  different  carriers 
and  their  associated  tones.  Care  must  be  taken  that  each  carrier  with  its 
tones  is  allowed  ±  5*  kHz  deviation  for  doppler. 

In  this  study  we  do  not  pursue  further  the  FDMA  option  since  it  does  not 
seem  amenable  to  making  the  receiver  cheap.  On  the  other  hand,  TDMA  has  some 
basic  advantages  as  discussed  in  Section  7. 

The  design  of  a  TDMA  Tone  Ranging  prototype  GPS  system  is  optimized  for 

the  following  general  parameters: 

The  carrier  frequency  is  -  1575  MHz. 

User  dynamics  are  generally  limited  to  one  g.  One  g 
dynamics  will  be  assumed  when  optimizing  system  para¬ 
meters  such  as,  for  example,  loop  bandwidth. 

Ranging  and  data  collecting  capability  must  be  present 
virtually  simultaneously,  i.e.,  both  functions  are  exercised 
during  the  dwell-time  on  a  single  same  satellite. 
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Basic  data  rates  and  basic  ranging  accuracies  are  the 
same  as  for  the  C/A  GPS  mode. 

Track  update  rates  must  be  higher  than  once  per  ten  seconds 
but  lower  than  once  per  second.  Indeed,  from  practical  ex¬ 
perience  with  tracking  civilian  aircraft,  update  intervals 
of  over  ten  seconds  seem  excessive  and  intervals  smaller 
than  one  second  do  not  greatly  contribute  to  better  performance. 

We  will  aim  generally  at  a  2  to  4  second  update  interval. 

In  the  remainder  of  this  section  we  give  a  synoptic  description  of  the 
prototype  GPS  tone  ranging  system,  followed  by  a  rationalization  of  some  im¬ 
portant  parameter  choices.  Although  it  is  concluded  that  a  very  adequate  GPS 
tone  ranging  system  could  be  designed,  it  is  not  concluded  that  it  promises 
great  simplicity  for  the  receiver  unless  real  breakthroughs  occur  in  reducing 
the  cost  of  phase  locked  loops. 

5.1  Tone  Ranging  System  Description 

In  this  section  we  summarize  the  salient  features  of  the  GPS  tone  ranging 
system  design:  the  power  budget,  the  chosen  tones,  the  expected  ranging  per¬ 
formance,  the  particular  choice  of  waveform. 

Table  5-1  shows  the  power  budget  for  the  GPS  tone  ranging  system.  The 
average  power  roughly  equals  that  of  the  GPS  C/A  mode.  The  TDMA  format  with 
twelve  time  slots  per  frame  is  used,  and  a  peak  radiated  power  twelve  times 
that  of  the  C/A  waveform  is  postulated.  This  is  5  dB  less  than  the  pulse  sys¬ 
tem  peak  power. 

Table  5-2  gives  the  clustered  ranging  tones.  Table  5-3  summarizes  the 
expected  range  error  budget.  The  estimated  ranging  error,  due  to  noise  only, 
derived  from  the  tones  is  28  ft.  The  resolvable  range  ambiguity  is  39.6  nmi. 

The  total  ranging  error  depends  largely  on  that  error  and  atmospheric  delays 
and  multipath.  The  latter  is  very  situation-dependent  and  has  not  been  studied 
further.  To  obtain  a  crude  estimate  of  user  position  error,  the  ranging  error 
must  be  multiplied  by  the  GDOP.  In  practice  a  factor  larger  than  GDOP  should  be 
applied  since  the  range  measurements  are  non-simultaneous.  The  tracker  (for 
example,  an  adaptive  Kalman-type  filter)  may  on  the  other  hand  reduce  the  result¬ 
ing  error  estimate  again  by  a  factor  typically  of  two  or  three,  depending  on  the 
user  motion  and  on  the  validity  of  the  tracker  equations  modeling  of  the  real 
world. 


The  principal  part  of  the  design  is  the  selection  of  the  particular  wave 
form.  It  seemed  necessary  to  adopt  a  waveform  at  least  similar  to  the  one 
shown  in  Fig.  5-1  where,  at  the  beginning  of  each  slot,  the  pure  carrier  is 
transmitted,  followed  by  a  time  interval  where  tones  are  transmitted,  then 
followed  by  an  interval  of  data  transmission.  The  remainder  of  the  timeslot 
is  a  guard  time  interval  protecting  against  overlapping  transmissions  of 
satellites  at  different  ranges.  The  advantage  of  a  short  "pure  carrier" 
interval  is  that  a  variable  bandwidth  PLL  can  quickly  acquire  and  lock  onto 
the  carrier.  From  the  doppler  on  the  carrier  the  precise  frequencies  of  the 
several  tones  can  then  be  predicted  so  that  lock-on  for  the  tone  PLL  at  base¬ 
band  can  occur  very  fast.  The  length  of  the  time  interval  for  tones  is  then 
mostly  determined  by  whether  the  tones  are  acquired  in  series  or  in  parallel. 
For  the  design  outlined  here  the  probability  of  bit-error  is  less  than  1.1  10 

Figure  5—2  shows  a  block  diagram  of  a  GPS  tone  ranging  receiver  the  para¬ 
meters  of  which  are  discussed  in  Section  5-2  and  Appendix  C. 


TABLE  5-1.  GPS  TONE  RANGING  POWER  BUDGET 


J  1  1  '  ”  —  1 

Radiated  Power  (25  w  average; 

12  x  25  w  in  TDMA)  300  w  peak 

25.0  dBw 

Satellite  Antenna  Gain 

11.5  dB 

Path  Loss  at  L^  ■  1575  MHz; 
distance  28524  km. 

-184.5  dB 

Atmospheric  Loss 

-1.0  dB 

Receiver  Losses  (incl. 
polarization  loss) 

-2.0  dB 

Receiver  Antenna  Gain 

0.0  dB 

Received  Power 

-151.0  dBw 

Thermal  Noise  (kT) 

-204.0  dB/Hz 

Noise  Figure 

5.0  dB 

-199.0  dBw/Hz 

Received  C/N 

0 

48.0  dBHz 

Margin 

3.0  dB 

Design  Value  for  C/Nq 

45.0  dBHz 

45 
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TABLE  5-3.  ERROR  BUDGET  OF  GPS  TONE  RANGING  SYSTEM 


Basic  Ranging  Error  cr^ 

28  ft 

Satellite  Ephemeris 

5  ft 

Atmospheric  Delay 

5-50  ft 

RSS  29-58  ft 

Multipath 


(variable) 


C:  20  msec,  CARRIER  ONLY 
T:  20  msec,  CARRIER  FM  MODULATED  BY  TONES 
D:  109.7  msec,  CARRIER  DPSK  MODULATED  BY  DATA 
G:  17  msec,  GUARD  TIME 


FIG.  5-1  TONE  RANGING  GPS  SIGNAL  FORMAT 
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FRONT  END  AMPLIFIERS 
MIXERS 
FILTERS 


FIG.  5-2  BLOCK  DIAGRAM  OF  GPS  TONE  RANGING  RECEIVER 
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5.2  Analysis  Defining  the  Main  System  Parameters 

'The  main  aspects  of  the  system  design  are  (1)  the  choice  of  the  waveform, 
i.e.,  the  determination  of  the  several  time  intervals  shown  in  Fig.  5.1;  (2) 
choice  of  the  tones  and  (3)  choice  of  RF  power  in  the  tones  or  equivalently  the 
choice  of  FM  modulation  indices  for  the  tones. 

5.2.1  Signal  Format 

A  signal  format  well  suited  to  the  multitone  ranging  problem  consists  of 
a  period  of  pure  carrier,  followed  by  a  period  of  carrier  FM  modulated  by 
the  several  sidetones,  followed  by  a  long  period  of  carrier  DPSK  modulated  by 
binary  data.  (See  Fig.  5.1).  The  length  of  the  various  time  intervals  was 
chosen  to  allow: 

acquisition  from  absolute  cold  start,  i.e.,  when  no  a  priori 
data  on  satellite  whereabouts  are  available  to  the  user; 

re-acquisition  within  the  time  interval  of  pure  carrier  after 
every  eleven  timeslots  of  silence,  coping  with  frequency  drift 
caused  by  relative  user-satellite  motion  and  oscillator  insta¬ 
bilities; 

acquisition  of  the  several  tones  and  settling  of  the  tone  PLL's 
and  read-out  of  the  individual  phases; 

conveying  a  sufficient  data  rate. 

This  section  discusses  further  the  choice  of  basic  cycle  time  T£  and  the 
length  of  time  blocks  C,  T,  D  and  G  of  Fig.  5.1. 

The  cycle  time  T,  determines  the  rate  at  which  the  user  position  is  up¬ 
dated  by  a  range  measurement  from  the  same  satellite.  We  want  T  to  be  on  the 

order  of  two  to  four  seconds.  Also,  T  determines  T  ,  the  timeslot  for  each 

c  s 

satellite  since  Tc  *  12  Tg.  The  guardtime  G  is  easily  determined  to  be  17  msec, 
the  time  difference  between  the  range  to  a  satellite  at  5  elevation  and  at  the 
zenith. 

5. 2. 1.1  Determination  of  the  Time  Interval  for  Radiating  Pure  Carrier 

The  time  interval  C,  during  which  pure  carrier  is  radiated  with  a  received 

C/N  of  45  dBHz,  must  be  long  enough  to  allow  phase-lock  to  occur.  The  frequency 
o 

excursion  which  the  PLL  has  to  cope  with  depends  on  the  length  of  time  the  phase 
tracking  for  the  particular  satellite  was  interrupted,  the  assumed  oscillator 
drifts  and  the  assumed  user  dynamics. 
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The  PLL  used  will  initially  work  with  a  large  bandwidth  BL*  The  pull-in 
frequency  threshold  Afp,  for  a  high  gain  second  order  loop  is  given  by 


For  a  short  enough  cycle  time  such  as  T  -  2  sec,  for  user  accelerations  not 

c 

more  than  a  few  g  and  for  an  oscillator  with  stability  10  Hz/Hz,  only  a  few 
hundred  Hz  of  frequency  offset  will  have  accumulated. 

For  Bt  ■  500  Hz,  the  pull-in  range  for  a  critically  damped  loop  is  roughly 

id 

±  400  Hz,  in  which  case  pull-in  occurs  within  5.5  msec.  If  the  frequency  off¬ 
set  is  less  than  225  Hz,  lock-on  occurs  in  about  1  msec. 


The  dynamic  lag  error  for  such  large  bandwidth  is  negligible  and  the 
statistical  phase  error,  E,  equals  about  8°  for  C/Nq  of  45  dBHz .  This  phase- 
angle  is  important  because  cos(E)  will  scale  the  power  available  in  the  PLL  of 
the  individual  tones  at  baseband.  For  this  reason  a  variable  bandwidth  PLL 
could  be  used,  with  small  bandwidth  being  used  once  lock-on  has  occurred. 
Choosing  time  interval  C  to  be  20  msec  for  a  cycle  time  Tc  of  2  seconds  there¬ 
fore  is  adequate. 

5. 2. 1.2  Determination  of  Time  Interval  T  for  Radiating  Tones 

In  tone  acquisition,  two  issues  are  of  concern.  The  duration  of  the  tones 
should  be  long  enough  for  the  individual  PLL's  to  settle  down.  The  second 
issue  is  the  final  phase  accuracy  with  which  the  phase  is  being  tracked. 

The  settling  Lime  again  depends  on  the  initial  phase  uncertainty  and 
on  the  noise  bandwidth  of  the  tone- PLL. 
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Since  we  have  a  fixed  set  of  tones,  and  the  carrier  on  which  the  tones 

are  modulated  is  phase-locked  before  the  tones  come  on,  the  frequency  offset 

on  the  tones  at  baseband  is  7—  f_,  where  f  denotes  the  tone  frequency  at  base- 

t  c  u  in 

band,  fc  the  carrier  frequency  and  fQ  the  doppler  offset  of  the  carrier.  For 
the  highest  tone  and  maximum  doppler  of  5  kHz,  this  frequency  offset  equals 
1.6  Hz.  In  principle  one  could  read  fQ  and  correct  the  VCO  of  the  tone  PLL. 
However,  the  noise  bandwidth  of  the  tone  PLL  will  be  large  enough  to  take  care 
of  1.6  Hz  initial  frequency  offset. 

For  a  tone  PLL  bandwidth  of  70  Hz  (tailored  to  a  tone-to-noise-power 
density  ratio  of  35  dBHz) ,  the  phase  error  is  about  10°  for  1  g  dynamics.  The 
PLL  will  then  lock  in  about  7  msec  for  a  32  Hz  offset.  However,  as  argued 
above,  the  offset  is  usually  much  smaller. 


In  order  to  accumulate  enough  signal 

needs  roughly  -r —  seconds  or  another 

“L, tone 

signal  integration  happen  simultaneously, 
interval  T  a  total  of  20  msec. 


for  a  lock-indicator  to  work,  one 
14  msec.  However,  the  lock-on  and 
consequently  we  allow  for  the  time 


5.2. 1.3  Choosing  Time  Interval  for  Transmitting  Data 

A  quick  calculation  shows  that  the  data- interval  will  have  to  be  consid¬ 
erably  longer  than  the  tone-interval.  Also,  to  keep  error  rates  down,  high 
C/Nq  is  desirable  for  a  given  bit  rate.  All  this  leads  us  to  have  a  separate 
data  interval  D,  and  to  DPSK  modulate  the  carrier  directly  during  that  Interval. 
Data  demodulation  circuitry  consists,  for  example,  of  a  Costas  loop  with  data 

detection.  For  C/N  ■  45  dBHz  and  about  11  dB  SNR  (for  bit  error  rate  10  ^ 
o 

using  DPSK  data  modulation),  a  data  rate  of  2.5  kbits  per  second  could  be 
supported. 

If  T£  is  chosen  to  be  2  sec,  the  number  of  bits  ,  R,  which  can  be  transmitted 
during  the  interval  ,  D,  where 
T 

D*jj  -  G  -  T  -  C  -  109.78  msec, 

is  137  bits.  This  is  sufficient  since  it  exceeds  the  50  bits  per  second  for 
the  present  GPS  design. 
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5.2.2  Choosing  Number  of  Tones  and  the  Tone  Frequencies 


"  The  number  of  tones  and  the  choice  of  frequencies  is  based  on  the  assumed 
initial  uncertainty  in  the  user's  position  (lowest  tone),  the  required  ranging 
precision  (highest  tone)  and  the  capability  of  the  tone  phase  lock  loops  to 
track  the  tone  with  a  given  phase  error  under  assumptions  of  moderate  user 
dynamics  and  of  available  input  signal  to  noise  power  for  the  tone  loops. 

The  tones  actually  modulating  the  carrier  are  the  sums  and  the  differences  of  the 
basic  tones  on  which  the  range  calculation  is  based,  and  are  "clustered"  in 
the  same  part  of  the  RF  spectrum.  FM  or  PM  modulation  of  the  carrier  by  the 
ranging  tones  seems  generally  preferable  to  AM  from  the  point  of  view  of  the 
transmitter.  The  RF  spectrum  in  that  case  is  generally  wide,  many  times  the 
baseband  of  the  modulating  tones.  In  a  tone  ranging  system  no  attempt  is  made 
to  "demodulate"  the  signal  and  collect  the  wideband  energy  efficiently.  In¬ 
stead  the  RF  spectrum  is  heterodyned  to  baseband  and  a  narrowband  filter 
(in  the  form  of  a  PLL)  is  applied  to  the  fundamental  spectral  line  of  each  tone. 
The  input  SNR  at  the  tone-PLL  is  a  function  of  how  the  tone-modulating  power 
is  distributed  among  carrier  and  the  tones.  The  distribution  is  determined 
by  the  modulation  indices.  Equality  of  indices  for  each  of  the  clustered 
tones  is  appropriate  in  our  design.  The  index  then  is  chosen  to  maximize 
the  power  in  the  fundamental  spectral  line,  with  as  little  power  as  possible 
in  intermodulation  spectral  lines  and  harmonics.  At  the  same  time  the  ratio 
of  power  in  the  fundamental  line  versus  the  total  power  is  optimized. 

In  the  following  paragraphs  some  rationale  is  given  for  choosing  the 
frequencies,  the  optimum  number  (four)  of  tones  and  the  optimum  index  of  70%. 

Tone  ranging  is  based  on  the  basic  principle  that  delay  in  propagation 
of  a  CW  signal  is  equivalent  to  phase  shift,  with  an  ambiguity  equal  to  the 
signal  period  (assuming  phase  can  be  measured  in  the  interval  [0,2ir]). 

If  the  a  priori  uncertainty  on  the  range  (due  to  imprecise  data  on 
satellite  position  and  user  position)  is  M  km,  then  the  lowest  frequency  tone 
must  have  a  frequency  no  higher  than  H  ’  where  c  is  the  velocity  of  light. 


By  tracking  the  tone,  the  phase  is  determined.  Let  sin(wLtg)  be  the  tone 
as  transmitted  by  the  satellite  (forgetting  the  carrier  modulation  and  demodu¬ 
lation  process).  Let  sinfc^  tR  +  <p)  be  the  received  tone,  with  phase  $.  Sym¬ 
bols  t  and  t_  are  time-readings  according  to  the  satellite  clock  and  the  re- 

S  In. 

ceiver  clock  respectively.  The  unambiguous  range  portion  is 


^L  [u>  (t  -t  )  +  <fi)  -  c(t_-t  )  +  \  •  <J>  . 
2tt  L  R  s  R  S  _ 


(3) 


The  portion  c(tR-tg)  will  be  provided  by  the  solution  of  four  ranges  to  four 
satellites  resulting  in  three  dimensional  user  position  and  clock  bias.  For 
simplicity  we  assume  the  satellite  times  of  the  four  satellites  are  perfectly 
synchronized  and  equal  to  true  time.  The  term  c(t_-t  )  represents  the  user 
clock  bias.  From  formula  (3)  we  see  how  the  range  error  is  related  to  the 
phase  error  and  the  wavelength: 


aR  “  V  2tt 


(4) 


In  principle,  if  the  phase  error  were  arbitrarily  small  the  range  error  would 
be  arbitrarily  small.  However,  a ^  is  limited  by  the  signal- to-noise  ratio  in  the 
PLL  and  the  presence  of  dynamics.  Therefore  higher  frequency  tones  must  be  chosen 
to  resolve  the  basic  uncertainty  area.  Finally  the  highest  frequency  tone  provides 
the  basic  ranging  accuracy  in  accordance  with  the  same  formula  (4) . 

We  continue  our  analysis  for  a  clustered  set  of  tones,  since  we  have  no 
intention  of  using  these  basic  frequencies  (f^..,,  f^) ,  but  want  to  use  clustered 
tones  (f'^,....  f"N)  instead  obtained,  for  example,  as  follows: 


f'l  *  fl 


f'2 


frf2 

frf3 


frfN 
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or  as  follows: 


V 

n 


f.-f 

1  n 


f.+f 
1  n 


for  n  even 
for  n  odd. 


For  tone  ranging  with  clustered  tones,  the  finest  range-reading  will  be 
improved  by  a  factor  or 


a 


R 


1  X 

7n”  <J>*  2tt  , 


(5) 


for  N  clustered  high  frequency  tones  with  average  wavelength  X,  when  each  tone  is 
tracked  with  a  phase  error  having  standard  deviation 


On  the  other  hand,  for  ambiguity  resolution,  if  the  tones  f'  are  tracked  with 
error  typified  by  the  corresponding  error  on  f^  is  /2  ,  or  worse  by  a 


factor  of  /2 


Table  5.4  shows  values  for  aD  according  to  formula  (5)  assuming  a  standard 
deviation  of  phase  error  of  10  degrees,  which  is  obtainable  with  a  PLL  with  in¬ 
put  SNR  of  35  dBHz  and  subject  to  1  g  dynamics.  The  SNR  figure  corresponds  to 
a  modulation  index  of  70%  and  is  justified  in  Appendix  B  on  Power  Distribution. 
For  aR  equal  to  30  ft  (the  most  interesting  design  value)  the  highest  frequency, 
which  is  also  roughly  the  average  cluster  frequency,  is  then  from  533  to  410  kHz, 
depending  on  the  number  of  tones  used. 


The  number  of  tones  is  determined  by  the  desired  probability  of  correctly 
resolving  all  range  ambiguities  from  lowest  to  highest  reconstituted  tones,  or 
equivalently  by  the  probability  that  any  of  the  N-l  ambiguities  was  falsely 
resolved.  The  latter  probability,  with  a  threshold  of  10  was  used. 


10  then  the  phase  error  on  the 


If  each  tone  f"  can  be  tracked  with  a, 

n  o  ** 

reconstituted  tone  f  is  /2~  a,  -  14  with  corresponding  ranging  error 


R,n 


2i 
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TABLE  5-4.  RANGE  ACCURACY  AND  HIGHEST  TONE 


°R  <ft) 

N-3 

fL(Hz) 

N=4 

f^Hz) 

N**5 

fj^CHz) 

10 

1.6  MHz 

1.4  MHz 

1.2  MHz 

30 

533  kHz 

463  kHz 

410  kHz 

100 

160  kHz 

140  kHz 

120  kHz 

300 

53.3  kHz 

46.3  kHz 

41  kHz 

*» 


If 


n-1 


k  then  the  probability  that 
n 


0D  >  X  . 
R,n  n-1 


resulting  in  incorrect  ambiguity  resolution  is 


,  /  Zir 
erf  (  — 


/2  a  ,k 
1  <p  n 


For  N  tones,  there  are  N-1  ambiguity  resolving  measurements;  and  the  overall 
probability  of  false  ambiguity  resolution  is 


p  -  1-  (1-P  )  (1-P  )  ...(1-P  ) 

FN  FN-1  F2  1 

For  equal  ratio  k  *  k  .  . . .  ■  a ,  one  has  P_  *  ?_  *  . . . ,  and  then 

fn  fn-i 

Pv  *  l-(l-p_  )N_1 
F  fn 

PF 

or  P  =  —L  for  small  Pp  . 


Table  5.5  shows  how  Pp  and  Pp  vary  as  a  function  of  N  and  the  assumed  for 
t'  of  10°.  Clearly  3  tones  are  too  few,  but  4  tones  are  satisfactory.  Although 
5  tones  may  seem  better  than  4,  it  constitutes  an  unnecessary  additional  com¬ 
plexity  (one  more  PLL  in  a  parallel  setup) .  Therefore  we  chose  the  number  of 
tones  to  be  equal  to  4. 

Reverting  to  our  previous  table  we  see  that  »  463  kHz  is  then 
the  indicated  highest  frequency.  The  lowest  frequency  f^,  as  pointed  out,  is 
determined  by  initial  range  uncertainty.  If  almanac  data  are  available  on  all 
satellites,  the  initial  range  uncertainty  is  dominated  by  the  user  position 
uncertainty,  say,  55  km.  Then  f^  -  6  kHz.  Having  determined  the  highest  and 
lowest  frequency  and  knowing  the  number  of  tones  we  determine  the  ratio  a: 


f 

f 


1 

4 


463 

6 


or  a  ■  4 , 26 
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The  tones  therefore  are 

f1  -  463  kHz,  f2  -  108.69  kHz,  f3  -  25.51  kHz,  f4  -  6  kHz. 

However,  we  now  adjust  this  set  of  frequencies  using  some  engineering  judgment. 
We  adjust  f^  upward  to  511.3  kHz  which  is  one  tenth  of  the  basic  GPS  oscillator 
frequency.  Since  all  tones  must  have  definite  phase  relationships  (for  example 
all  phases  are  zero  at  t  ■  0)  it  is  easier  to  have  simple  tone  frequency  ratios 
For  example,  a  ■  5  gives  “  102.26,  f^  *  20.452,  f^  *  4.0904  kHz  with  finest 
accuracy  27.3  ft.  (8.2  m)  and  resolvable  initial  uncertainty  of  39.6  nmi 
(73.4  km)  and  probability  of  falsely  resolving  ambiguity  of  1.1  x  10 

The  tones  actually  modulating  the  carrier  are  then,  for  example: 

£'  -  511.3  kHz 

f'  -  409.04  kHz 

f'  -  490.85  kHz 
t'  -  507.21  kHz  . 

Finally  we  discuss  the  choice  of  modulation  index.  In  Fig.  5-3  it  is  shown 
that  the  suppression  factor,  which  is  defined  as  the  ratio  of  power  in  the  tone 
spectral  line  to  the  total  carrier  power,  has  a  maximum  for  the  FM  modulation 
index  A  equal  to  .7.  The  SNR  at  the  input  of  the  tone  PLL,  for  C/Nq  *  45  dBHz, 
is  then  about  38  dBHz.  The  SNR  available  for  the  tone  PLL  actually  used  in 
our  design  was  35  dBHz,  allowing  for  a  3  dB  margin. 


r  K  (08)  SUPPRESSION  FACTOR 


FIG.  5-3  SUPPRESSION  FACTOR  AS  FUNCTION  OF  FM  MODULATION  INDEX 
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5.3  Discussion  on  Tone  Ranging  Performance 

%The  main  conclusion  on  cone-ranging  as  an  alternate  waveform  is  that  one 
could  indeed  design  an  efficient  tone-ranging  system,  as  we  have  demonstrated. 

Multitone  systems  without  PN  encoding  are  susceptible  to  multipath  effects. 
It  is  possible  however,  that  for  a  clustered  tone  system,  multipath  effects 
could  be  averaged  out.  This  issue  requires  further  study. 

A  multitone  system  without  PM  encoding  could  be  easily  jammed,  since  all 
that  is  necessary  to  disrupt  it  is  to  reproduce  one  or  several  tones. 

The  multitone  system  designed  here,  with  periods  of  pure  carrier  being 
radiated,  with  some  additional  acquisition  circuits,  could  "cold-start"  more 
easily  then  the  present  C/A  GPS  system.  A  cold-start  condition  corresponds  to 
total  lack  of  a  priori  information  on  the  satellite  positions.  The  only 
knowledge  assumed  is  a  crude  estimate  of  user  position  (error  less  than  40 
nmi) .  Appendix  B  provides  some  details  on  a  proposed  method  for  cold-starting 
the  multitone  ranging  system. 

Tone  ranging  was  found  to  be  practical  only  in  a  TDMA  context.  Its 
ranging  and  data  carrying  capability  can  be  favorably  compared  to  those  proper¬ 
ties  of  the  GPS  C/A  system.  As  an  alternate  waveform  to  the  GPS  C/A  waveform 
it  is  generally  more  complex  and  is  inferior  to  the  pulsed  system  of  Section  4, 
but  its  peak  power  requirements  are  lower. 

Finally,  we  present  Table  5-6  with  a  summary  of  the  main  features  and 
performance  characteristics  of  Tone-Ranging  as  presented  by  our  chosen  design. 

The  final  criterion,  cost,  is  not  included.  Based  on  our  conceptual 
design  of  a  multitone  ranging  system,  there  is  little  indication  that  it  would 
cost  less  than  a  sequential  GPS  set  due  to  the  large  number  of  phase- locked 


Power  Levels 

Transmitted  Average  Power 

Transmitted  Peak  Power 

Transmission 

Power  at  Receiving  Antenna 

Received  C/Nq(5  dB  noise  figure) 

25W 

25  dBW 

TDMA  format  150  msec 
burst  per  2  sec. 

-151  dBW 

48  dBHz 

Signal  Design  Features 

Ranging  Capability  Provided  by 

Set  of  4  tones 

Data  Carrying  Method 

DPSK  of  Pure  Carrier 

Data  Rate 

Burst  of  137  bits  in 

109  msec  or  66  bps  avg. 

Data  Repetition  Period  for  User 

24  sec 

Receiver  Features 

Ranging  Loop 

Bank  of  4  PLL's 

Carrier  Tracking  Loop 

Costas  PLL 

Data  Demodulation 

Coherent 

Start-Up  Circuitry 

Yes 

Cold  Start  Capability 

Easy 

Range  Update  Rate 

0.5  Hz 

Measurement  Type 

R,R 

Sensitivity  to  Multipath 

Sensitive 
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6 .  PRFDMA  WAVEFORM 


In  the  class  of  pulsed  waveforms,  where  the  short  high  energy  pulse  is 
adopted  as  the  basic  building  block,  there  is  a  large  variety  of  useful  signal 
structures.  One  such  was  the  TDMA  Pulsed  Design  of  Section  TV.  The  different 
parts  of  the  waveform  were  designed  to  correspond  to  specific  functions:  a 
pulse  burst  for  ranging  and  synchronization  followed  by  PPM  for  data  trans¬ 
mission  at  a  much  lower  duty  cycle. 

Although  a  short-term  duty  cycle  is  a  good  choice  for  a  simple  receiver, 
it  puts  demands  on  the  transmitter.  Another  pulsed  waveform,  which  relaxes  the 
requirements  on  the  transmitter,  consists  of  a  simple  pulse  train  at  a  fixed 
PRF.  Such  a  waveform  will  also  reduce  the  demands  on  the  user's  clock  in  terms 
of  stability.  Each  antipodal  satellite  pair  uses  a  slightly  different  PRF  to 
avoid  repeated  pulse  overlap  at  the  receiver,  hence  the  name  applied  to  this 
waveform. 

A  quick  design  study  of  one  such  PRFDMA  waveform  was  made  and  is  summarized 
below.  The  general  conclusion  was  that,  if  computation  costs  keep  decreasing 
at  the  rate  experienced  in  the  last  decade,  such  a  waveform  may  compete  favor¬ 
ably  with  the  TDMA  pulsed  system,  because  of  the  simplicity  of  the  RF  portions 
in  transmitter  and  receiver.  It  was  found  to  relax  clock  requirements  but  to 
be  inferior  in  acquisition  capability. 

6 . 1  Ranging 

The  PRFDMA  format  is  a  signal  format  in  which  all  satellites  transmit  a 
single  pulse  at  a  given  PRF,  rather  than  a  burst  of  pulses  over,  say,  20  msec, 
repeated  every  .5  seconds  as  described  in  Section  4  .  Figure  6-1  shows 
a  PRF  waveform  choice.  Assuming  a  2  kW  satellite  transmitted  peak  power,  the 
received  power  level,  is  56.1  dBHz.  For  the  chosen  pulse  length  of  41  usee, 
this  corresponds  to  a  SNR  of  12.2  dB  when  no  frequency  offset  is  present.  A 
5  kHz  frequency  offset  causes  the  SNR  to  drop  to  11.6  dB  and  a  9.2  kHz  offset 
results  in  a  SNR  of  10  dB.  This  is  the  minimum  detection  threshold  for 
operation. 
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164  CHIPS 


FIG.  6 


PRF.  =  240+i  FOR  1.6  SEC 

=  240  -  i  FOR  NEXT  1.5  SEC 


1  SKETCH  OF  PRFDMA  WAVEFORM 


The  pulse  length  of  41  ysec  was  chosen  to  allow  tolerance  to  a  frequency 
offset  of  up  to  9.2  kHz  from  all  sources  (relative  user-satellite  motion  and 
clock  frequency  offset  or  drift),  while  still  guaranteeing  a  SNR  of  10  dB  for 
pulse  detection.  This  pulse  is  band-spread  by  means  of  a  4  MHz  BPSK  modulation. 


The  ranging  accuracy  is  derived  from  the  bandspreading,  using  PN  encoding 
with  a  164-chip  code,  identical  for  all  satellites.  Ranging  based  on  a  CCD 
matched  filter  implementation  is  expected  to  result  in  a  standard  deviation 
on  a  single  pulse  measurement  of  ap  -  ,  where  A  is  the  chip  length  of 

259  nsec  or  -  75  ft 


m 


Pseudorange  accuracy  is  not  based  on  a  single  pulse, 
but  on  the  output  of  a  pseudorange  tracker  updated  by  range  readings  at  the 
PRF  rate. 


A  ranging  pulse  is  received  approximately  every  4.2  msec.  The  time  of 
transmission  of  these  pulses  (based  on  the  satellite  ID  and  the  data)  is  also 
known.  (This  may  involve  counting  the  number  of  pulses  with  respect  to  a 
reference  point  in  time  such  as  the  Z-count.)  One  can  then  track,  for  example 
with  a  simple  fixed  gain  (a, 8)  tracker,  the  received  pulse  time  t  versus 
satellite  transmit  time  s.  The  tracking  allows  very  accurate  prediction  of 
the  arrival  time  of  the  next  pulse,  thus  helping  in  associating  the  correct 
pulse  with  the  correct  track,  controlling  the  false  alarm  rate.  But  more  im¬ 
portantly  In  the  context  of  accuracy,  tracking  increases  ranging  accuracy  and 
allows  strobing  of  all  pseudorange  tracks  at  a  common  time  t  on  the  user's 
clock  and  determination  of  corresponding  time  s  at  which  each  satellite  would 
have  sent  a  pulse  arriving  at  the  user  at  time  t.  Correct  user  position  deter¬ 
mination  can  then  occur,  corresponding  to  time  t  even  though  the  readings  from 
the  several  satellites  are  asynchronously  received. 

The  twelve  PRFs  needed  could  be  defined  by  (PRF^  ■  240  ±  i,  for  i  -  0,11) 
where  the  minus  value  holds  for  1.5  seconds  and  the  plus  value  over  the  next 
1.5  seconds  for  an  overall  average  of  240  pulses  per  second  over  three  seconds. 
Obviously  with  as  many  as  eleven  satellites  visible,  pulses  will  occasionally 
overlap  (about  one  out  of  ten) .  However  they  must  overlap  to  within  two  chips 
for  them  to  overlap  undistinguishably  (probability  of  about  one  out  of  eight 
hundred).  Such  a  condition  can  be  detected  and  no  decision  made  on  data  or 
ranging  information  derived  from  that  single  pulse. 
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It  is  easy  to  show  thac  for  che  above  PRF  races  two  indistinguishably 
overlapping  pulses  (aligned  within  two  chips)  from  two  satellites  will,  for 
the  next  pulse  interval,  be  separated  by  about  68  chips.  For  a  signal  which 
had  already  been  acquired,  the  pseudorange  tracker  would  simply  coast  for  one 
update.  The  data  demodulation  would  be  based  on  the  remaining  "clear"  pulses. 

The  pseudorange  tracker  helps  to  sort  the  pulses  into  the  right  satellite 
file,  the  identity  of  which  is  re-affirmed  once  per  half  second  by  the  data. 

If  any  contradiction  in  observed  PRF  and  the  identification  exists,  part  of 
che  data  may  have  to  be  purged  or  re-acquisition  procedures  initiated. 

6.2  Data  Demodulation 

The  pulses  are  PN-encoded.  One  may  choose  a  different  code  to  represent  a 
zero  or  a  one  or  use  PPM  for  data  modulation.  Pulse  position  data  modulation 
would  increase  che  computation  load  and  make  initial  startup  more  complicated. 

Code  division  increases  the  hardware  complexity.  One  bit  would  be  transmitted 
per  four  pulses  for  a  data  rate  of  60  bps  which  corresponds  to  the  50  bps  of 
GPS,  and  another  10  bits  of  additional  data,  for  example,  for  extra  identification. 

6.3  Satellite  Identification 

The  PRF's  of  che  several  satellites  are  different;  but  in  order  to  main¬ 
tain  synchronism  with  the  P-signal  at  regular  intervals  (for  handover  and  for 
simultaneous  change-over  to  refreshed  data  once  every  hour)  each  satellite  can 
be  assigned  two  PRF's,  averaging  to  the  same  value  for  all  satellites,  as  was 
suggested  above  in  Section  6.1. 

6.4  Clock  Stability  and  Computer  Load  Implications 

In  Section  9  an  analysis  is  presented  on  the  implications  of  the  PRF 
update  rate  of  4  msec  for  this  approach  versus  .5  second  for  the  pulsed  and 
tone-ranging  TDMA  candidates,  on  the  clock  stability  requirements. 

With  the  user-clock  timing  error  contribution  to  ranging  restricted  to 
13  nsec  and  for  the  assumption  of  white  noise  the  clock  stability  requirement 
expressed  in  Allan  variance  can  be  relaxed  roughlv  by  a  factor  of  five  when 
using  a  4  msec  update  interval  instead  of  .5  seconds.  In  terms  of  computer 
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load,  it  was  judged  necessary  to  have  dedicated  computer  chips,  estimated  to 
cost' a  few  dollars  each,  to  perform  pseudorange  tracking  in  order  to  free  the 
microprocessor  for  other  tasks. 

In  conclusion,  a  PRFDMA  system,  with  relaxed  requirements  on  the  satellite 
transmitter  and  the  user  clock,  can  be  designed  with  satisfactory  performance. 
The  receiver  would  be  simple  but  the  acquisition,  pulse  sorting  and  pseudorange 
tracking  tasks  put  additional  burdens  on  the  microprocessor  as  compared  with 
the  TDMA  pulsed  design. 
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7.  IMPLICATIONS  OF  THE  TDMA  FORMAT 


In  this  section  we  discuss  some  advantages  of  a  TDMA  format  for  sharing 
the  channel  capacity  amongst  the  24  GPS  satellites. 

The  low-cost  GPS  receiver  will  most  likely  be  sequential,  with  the  receiver 
cycling  between  a  minimum  of  four  signals  needed  for  position  and  clock  bias 
calculation,  and  possibly  tracking  a  fifth  or  sixth  signal  as  well,  to  cope 
with  signal  interruptions  caused  by  maneuvers  or  to  prepare  for  re-selection 
of  satellites  to  maintain  acceptable  geometric  dilution  of  precision. 

In  che  case  of  four  satellites,  only  during  252  of  the  time  that  signals 
are  continuously  transmitted,  are  they  actually  used  by  the  sequential  receiver. 
In  che  proposed  TDMA  format,  a  satellite  transmits  only  during  one-twelfth  of 
the  time,  in  a  fixed  timeslot.  For  the  same  average  transmitted  power  the 
received  signal  level  can  be  about  10  dB  stronger.  Also,  it  can  be  received 
during  1002  of  the  time  that  it  is  transmitted  and  since  satellites  transmit 
only  in  assigned  timeslots,  separated  by  a  time  interval  to  avoid  signal  over¬ 
lap.  the  satellite  signals  do  not  mutually  interfere. 

The  advantages  for  data  collection  are  obvious.  In  the  case  of  GPS  C/A 
and  a  sequential  receiver,  message  portions  are  received  piecemeal  and  out  of 
sequence.  In  addition  to  this,  it  is  hard  to  match  the  cycle  time  of  the  re¬ 
ceiver  against  a  particular  subdivision  of  the  messages.  Due  to  the  differen¬ 
tial  delay  between  satellites  (worth  up  to  1  bit)  and  re-acquisition  times 
(worth  several  bits) ,  several  data  bits  between  pieces  of  the  message  will 
most  probably  be  lost. 

In  the  proposed  TDMA  format  there  is  no  difficulty  in  assembling  the 
message  portions  in  sequence  and  it  is  done  without  loss  of  bits. 

The  counter-argument  could  be  that,  although  now  transmission  and  reception 
are  one  hundred  percent  matched,  the  receiver  is  idle  except  during  four  out  of 
twelve  timeslots,  or  752  of  the  time.  And  since  dat3  are  received  in  a  burst, 
the  higher  instantaneous  data  rates  will  imply  a  more  powerful  and  therefore 
more  costly  processor.  However,  the  latter  argument  does  not  hold,  since  the 


68 


circuitry  for  data  demodulation  for  the  pulsed  system  discussed  here  is  simple, 
and  unaffected  by  the  high  instantaneous  data  rate.  The  proposed  instantaneous 
rate  of  about  4  kbits  per  second  in  the  particular  TDMA  pulsed  design  of  Sec¬ 
tion  IV  is  within  the  capabilities  of  present-day  processors.  Let  us  point  out 
that  the  4  kbps  rate  is  a  consequence  of  our  choice  of  a  three  times  higher  average 
data  rate  and  should  be  correspondingly  reduced  when  comparing  the  receivers  for 
equal  performance. 

To  understand  better  the  discussion  on  GDOP  which  follows,  we  include  some 
general  conclusions  concerning  our  short  study  on  satellite  selection  based  on 
optimal  or  sub-optimal  GDOP. 

A  study  on  ephemeris  behavior  of  GPS  satellites  and,  in  particular,  of 
a  constellation  of  four  satellites  selected  for  their  good  GDOP  properties 
revealed  that: 

obtaining  the  "optimal"  selection  is  a  laborious  task 

optimality  is  usually  short-lived  (on  the  order  of 
ten  minutes) 

the  change  of  GDOP  as  a  function  of  time  for  a  selected 
constellation  can  be  quite  fast,  leading  possibly  to  a 
50%  deterioratation  in  a  time  span  of  ten  minutes. 

Often  the  best  constellation  is  one  that  is  made  up  of  three  low  elevation  satel¬ 
lites,  as  equally  distributed  as  possible,  together  with  one  high  elevation 
satellite.  Usually,  at  least  one  is  setting.  Good  GODP  and  longevity,  therefore, 
do  not  go  hand-in-hand. 

The  problem  is  that  when  a  single  satellite  becomes  undesired  (sets 
below  5°  EL  for  example) ,  it  is  most  often  true  that  no  other  satellite  can 
replace  it  to  maintain  a  satisfactory  GDOP.  Usually  a  completely  new  constel¬ 
lation  must  be  acquired  and  tracked.  The  short-livedness  (order  of  10  minutes) 
of  a  good  constellation  makes  this  a  difficult  problem  for  the  receiver.  The 
proposed  TDMA  system  would  not  have  this  problem. 
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For  the  pulsed  TDMA  system  (Section  4),  acquisition  and  cracking  of  a 
new  satellite  is  virtually  instantaneous  by  design.  The  signals  of  all  visible 
satellites  are  available  and  can  be  tracked  since  they  arrive  in  a  strict  non¬ 
overlapping  time-ordered  basis  at  the  receiver.  All  signals  have  virtually  equal 
strength  and  are  modulated  by  the  same  PN  code  and  are  on  the  same  frequency. 

There  is  no  adjusting  or  re-configuring  necessary  on  the  part  of  the  receiver. 

The  752  of  the  cycle  time  during  which  the  receiver  is  idle,  when  com¬ 
pared  to  the  C/A  GPS  receiver  tracking  only  four  satellites,  can  be  exploited 
to  great  operational  advantage. 

The  idea  of  using  four  range  measurements  to  four  selected  satellites 
based  on  their  associated  "good"  GDOP  can  be  totally  dispensed  with.  The 
algorithmic  complexity  of  selecting  and  continuously  reselecting  an  optimal 
or  at  least  sub-optimal  set  as  often  as  every  ten  minutes,  and  the  associated 
interruptions  for  re- acquisition  of  the  new  signals,  are  then  unnecessary. 

Also  the  problem  of  being  forced  into  unscheduled  re-selection  when  one 
satellite  signal  disappears  (sudden  turn)  need  never  occur. 

One  may  use  ranging  information  from  all  available  (visible)  satellites. 

When  more  than  four  ranges  are  used,  position  can  be  determined,  for  example, 
in  a  least  squares  sense  instead  of  a  deterministic  sense.  The  resulting 
position  coordinates  are  then  used  to  update  the  user  tracker.  In  many 
designs  it  is  natural  that  range  measurements  are  obtained  at  a  higher  rate 
than  is  necessary  for  the  user  tracker.  One  may  then  smooth  pseudorange 
measurements  before  combining  the  pseudoranges,  and  extrapolate  or  interpolate 
the  measurements  to  a  common  time  at  which  user  position  is  determined. 

It  should  be  pointed  out  that  while  the  GDOP  resulting  from  the  use  of 
all  visible  satellites  is  always  truly  optimal  (and  better  than  the  optimal 
GDOP  resulting  from  a  subset  of  only  four  visible  satellites) ,  no  such  claim 
can  be  made  in  general  for  the  GDOP  based  on  only  four.  Indeed,  it  is  almost 
unthinkable,  for  example,  that  with  eleven  visible  satellites  all  330  combinations 
of  four  satellities  would  be  selected  and  compared  based  on  their  GDOP.  Instead 
some  algorithm  will  first  reduce  the  number  of  combinations  to  a  tractable 
number  (maybe  ten)  and  only  then  compare  the  remaining  combinations  based  on 


70 


explicitly  calculated  GDOP  values.  Such  algorithms  are  therefore  at  best 
sub-optimal.  Also,  the  selected  combination  out  of  the  reduced  set,  based  on 
the  best  GDOP,  is  only  best  at  a  given  instant  of  time. 

It  is  quite  complicated  to  predict  how  long  the  selected  combination 

will  remain  best.  In  general,  one  could  compare  the  remaining  combinations 
based  on  the  expected  lifespan  of  the  combination  or  on  the  average  of  its 
GDOP  values  now  and  at  some  future  time.  This  makes  the  procedure  of  selection 
of  a  combination  even  more  sub-optimal. 

And  finally,  it  should  be  pointed  out  that  GDOP  is  a  scalar  measure 

expressing  linear  amplification  of  the  error  on  a  single  range  measurement, 

as  expressed  by  the  User  Equivalent  Ranging  Error  (UERE) .  In  reality  the 
geometric  dilution  is  more  pronounced  in  some  directions  than  in  others,  and 
these  directions  and  the  amplifications  suddenly  change  at  the  moment  of 
re-selection  of  one  or  more  satellites,  even  if  the  scalar  value  for  GDOP 
remains  the  same. 

If  all  visible  satellites  are  used,  the  GDOP  transitions,  in  scalar 
value  or  in  the  directionality  of  the  main  component  of  error,  will  be  less 
severe  since  satellites  are  replaced  automatically  one  at  a  time.  Although 
the  tracking  algorithms  are  exercised  more  often,  a  sizable  portion  of  soft¬ 
ware  associated  with  satellite  (re-) selection  and  the  time  to  exercise  it, 
can  be  saved. 
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8.  HANDOVER  FROM  ALTERNATE  WAVEFORM  MODE  TO  P  MODE 

vIt  is  possible  that  a  user  of  the  pulsed  C/A  waveform  may  wish  to  acquire 
the  normal  GPS  P-signal  for  increased  accuracy  of  ranging.  For  this  reason  the 
pulsed  signal  should  contain  handover  information  equivalent  to  that  provided 
in  the  normal  GPS  C/A  signal.  The  basic  requirement  is  the  inclusion  of  a 
data  word  containing  the- so-called  Z-count  and  a  means  of  recognizing  the  in¬ 
stant  of  time  at  which  this  Z-count  applies.  The  latter  part  of  this  require¬ 
ment  implies  a  basic  synchronism  between  the  P-waveform  and  the  pulsed  C/A 
signal,  as  discussed  in  more  detail  below. 

The  P-signal  is  formed  from  two  shift-register  sequences,  each  of  which  is 
short-cycled  to  15,345,000  chips.  Since  the  P-signal  chipping  rate  is  10.23  MHz, 
the  cycle  corresponds  to  precisely  1.5  seconds,  called  a  super-epoch.  At  GPS 
zero  time  (midnight  Saturday  night),  the  two  codes,  which  are  added  (modulo  2) 
to  form  the  P-code,  are  started  together  in  the  all-ones  state.  In  satellite 
Number  1  they  are  clocked  together,  while  in  satellite  Number  2  the  second  code 
is  delayed  one  chip,  and  so  on  up  to  a  delay  of  23  chips  for  the  last  satellite. 
At  the  end  of  the  first  super-epoch,  the  second  code  is  held  fixed  (clock 
pulses  are  inhibited)  for  32  chips  and  then  again  clocked  in  synchronism  with 
the  first  code.  The  second  code  marks  time  in  this  way  for  32  chips  at  the 
start  of  each  super-epoch,  thus  precessing  with  the  first  code  to  form  a  very 
long  sequence  for  the  resultant  P-code. 

One  week  contains  403,200  super-epochs,  hence  the  second  code  has  slipped 
by  12,902,400  chips  by  that  time,  after  which  the  second  code  is  reset  to  all- 
ones  and  the  entire  week-long  code  (a  different  one  for  each  satellite)  is 
repeated. 

The  Z-count  is  the  number  of  super-epochs  which  have  elapsed  since  zero¬ 
time,  and  it  requires  a  19-bit  word  for  full  representation.  The  normal  C/A 
signal  contains  data  (synchronous  with  the  C/A  PN  code)  organized  into  6-second 
subframes.  Each  subframe  begins  exactly  at  the  start  of  a  super-epoch  and  lasts 
for  four  super-epochs.  Each  subframe  also  contains  a  Z-count  word  of  17  bits 
(the  two  lowest  order  bits  are  truncated)  which  identifies  the  super-epoch  co¬ 
incident  with  the  start  of  the  next  forthcoming  subframe.  The  user  who  has 
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acquired  this  information  can  therefore  set  the  two  code  generators  which  will 
form'bis  P-signal  to  the  all-ones  state  and  know  exactly  how  many  chips  to  in¬ 
hibit  the  second  after  starting  the  first  at  the  time  of  the  next  subframe 
start.  He  will  be  in  step  with  the  super-epochs  of  the  arriving  P-code  to  an 
accuracy  of  a  fraction  of  a  C/A  code  chip,  i.e.,  to  an  accuracy  of  fewer  than 
ten  P-code  chips.  Each  of  the  codes  which  form  the  P-signal  can  be  slewed 
ahead  or  behind  by  a  method  analogous  to  the  slewing  of  these  codes  relative 
to  one  another  to  complete  the  P-code  acquisition. 

Any  C/A  waveform  can  provide  the  same  handover,  provided  it  is  synchronous 
with  the  P-code  super-epochs  in  some  way  and  also  provides  Z-count  Information. 
For  example,  in  the  pulsed  design  described  above  it  is  only  necessary  to  make 
the  0.5  second  frames  synchronous  with  the  P-code  super-epochs  so  that  every 
third  frame  begins  exactly  at  the  start  of  a  super-epoch.  A  19-bit  full  Z-count 
word  can  be  transmitted  during  the  triplet  of  frames  corresponding  to  one 
super-epoch.  Since  each  satellite  begins  transmitting  at  the  start  of  a 
precisely-timed  time  slot,  the  user  can  preset  his  P-code  generators  for  the 
start  of  the  next  forthcoming  super-epoch  and  also  allow  for  the  delay  expected 
between  this  super-epoch  start  and  the  arrival  of  the  sync  burst  from  the 
satellite  in  question  in  the  first  frame  of  this  super-epoch.  There  is  no 
requirement  for  a  precise  user  clock  to  span  a  long  interval  of  time  in  this 
procedure,  since  all  the  burden  for  precision  of  timing  (e.g.,  time-slot  start) 
is  placed  upon  the  satellite. 
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9.  CLOCK  STABILITY  REQUIREMENTS  IN  TDMA  SYSTEMS 


In  a  conventional  multichannel  GPS  receiver,  pseudorange  measurements 
from  all  tracked  satellites  are  continuously  available.  By  simultaneous 
sampling  of  these  outputs,  the  range  differences  necessary  for  user  position 
determination  can  be  obtained  without  reference  to  a  clock.  However,  in 
any  system  with  intermittent  transmissions  from  the  satellites,  such  as  the 
TDMA  systems  described  above,  pseudorange  for  each  satellite  is  obtained 
only  at  intervals,  and  must  be  timed  and  tracked  in  some  way  so  that  simul¬ 
taneous  measurements  can  be  synthesized  from  the  track  for  positioning. 

This  process  of  tracking  and  estimation  brings  with  it  a  dependence  on  clock 
behavior  and  a  requirement  for  short-term  clock  stability. 

In  the  PRFDMA  system,  pseudorange  tracking  is  necessary  to  sort  the 
pulses  from  different  satellites  and  also  to  perform  smoothing,  or  integration, 
of  the  pulses  from  each  satellite  in  order  to  build  up  the  accuracy  of  the 
individual  pseudorange  measurements.  A  simple  "alpha-beta"  tracker  has  been 
suggested  for  this  purpose,  one  for  each  satellite  in  track,  and  these  tracker 
outputs  can  be  used  to  predict  each  pseudorange  ahead  to  a  common  measurement 
time.  The  predicted  values  are  then  used  in  a  conventional  user  position 
calculation,  just  as  though  they  had  been  obtained  by  "strobing"  (simultan¬ 
eously  sampling)  a  set  of  continuous  pseudorange  outputs. 

In  the  TDMA  schemes,  the  situation  is  a  bit  different,  since  pseudorange 
is  obtained  once  per  frame  (one-half  second  interval)  from  each  satellite, 
and  every  measurement,  being  made  from  a  pulse  burst  or  tone-ranging  signal, 
already  has  the  accuracy  required  for  positioning.  The  measurements  are 
made  at  staggered  times,  by  the  basic  nature  of  the  TDMA  format,  and  require 
interpolation  (instead  of  smoothing  and  prediction)  to  a  common  measurement 
time  for  input  to  the  position  calculation  algorithm. 

In  all  cases,  the  pseudorange  measurements  actually  used  for  location 
are  estimates,  based  on  timed,  intermittent  data,  and  hence  carry  errors  due 
to  clock  imperfections  in  addition  to  the  usual  errors  of  TOA  (time-of- 
arrival)  measurement.  A  general  analysis  has  been  carried  out  of  the  effect 


of  clock  instability  on  user  position  accuracy  for  several  tracking  and  inter¬ 
polation  schemes,  and  for  two  standard  models  for  clock  behavior.  The  general 
conclusion  is  that  clock  timing  errors  are  less  serious  in  their  effect  on 
positions  than  TOA  errors  of  equal  magnitude,  due  to  the  correlated  nature  of 
clock  errors  and  the  tolerance  of  the  tracker-interpolation  schemes  to  low- 
frequency  clock  noise. 

In  this  section,  only  the  TDMA  case,  with  simple  two-point  interpolation, 
will  be  discussed  in  detail.  This  is  done  because  the  TDMA  systems,  with 
their  relatively  long  intervals  between  measurements ,  as  compared  to  the  PRFDMA 
system  place  the  greatest  demands  on  the  clock.  It  also  turns  out  that  the 
analysis  is  considerably  simpler  than  that  of  the  alpha-beta  tracker  model, 
for  which  results  will  be  stated  without  proof. 

Taking  the  pulsed  TDMA  system  first,  consider  the  ranging  bursts  from 
a  single  satellite,  and  let  the  measured  arrival  times  of  these  bursts  be 
R^,  according  to  the  user's  clock.  From  the  data  sent  by  this 

satellite,  the  user  knows  the  corresponding  transmission  times,  S^,  S2 . 

of  the  bursts,  according  to  the  satellite's  clock.  Both  times  are  adjusted 
to  refer  to  a  common  fiducial  mark  within  the  pulse  burst,  and  no  satellite 
clock  error  is  assumed.  In  effect,  the  user  can  plot  transmission  time, 

Sn,  versus  reception  time,  Rn>  and  join  the  points  with  straight  lines  for 
each  of  the  satellites  being  used.  When  position  computation  is  desired, 
a  single  time,  R,  i3  chosen  and  the  corresponding  transmission  times  are 
read  off  of  each  plot,  providing  simultaneous  pseudorange  for  location. 

Position  computation  is  asynchronous  with  the  TDMA  frame,  being  exercised  on 
demand  from  the  user  position  tracker. 

Since  the  data  points  were  joined  by  straight  lines,  the  process 
represents  simple  linear  interpolation,  for  each  satellite,  between  the  two 
measurements  which  straddle  the  measurement  time,  R.  Linear  interpolation 
is  appropriate,  since  smoothing  by  curve-fitting  is  not  required  to  improve 
accuracy,  and  the  effect  of  acceleration  in  user-to-satellite  range  (which 
is  dominated  by  the  user's  own  acceleration)  is  negligible  over  the  half- 


75 


second  intervals  involved  here,  amounting  to  no  more  than  one  foot  per  g  of 
acceleration. 

In  practice,  it  is  only  necessary  to  save  the  two  most  recent  measure¬ 
ments  (S  and  R  for  two  successive  frames)  for  each  satellite  in  ping-pong 
fashion  until  a  position  update  is  required.  If  this  request  comes  at  time 

R,  the  data  in  store  is  retrieved,  and  the  measurement  time  R  *  R  -  A  is 

o 

established,  where  A' is  the  frame  duration  for  the  system.  A  moment's  thought 

will  show  that  the  stored  data  will  contain  one  pair  with  reception  time 

earlier  than  R  and  one  pair  later  than  R  ,  for  each  satellite.  If  these 
o  o 

pairs  for  one  satellite  are  called  R^ ,  and  Rj,  S2  *  then  the  transmission 
time  estimate  corresponding  to  the  chosen  reception  time  Rq  is  simply 


A 

s 


o 
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It  might  seem  simpler  to  use  the  same  data  to  predict  ahead  to  time  R,  but 
the  general  analysis  shows  that  clock  error  effects  are  minimized  by  using 
RQ  in  the  middle  of  the  data  points,  and  that  prediction  outside  the  measure¬ 
ment  interval  causes  significant  degradation.  The  fact  that  the  position 
fix  is  older,  by  A  (.5  sec.),  than  otherwise  would  be  the  case,  is  of  no 
importance  since  the  user  position  tracker  will  be  predicting  position  in 
any  case  in  order  to  drive  the  user's  display. 

Suppose  that  true  time  is  denoted  by  t,  and  that  the  corresponding 
time  by  the  user's  clock  is 

R(t)  -  t  +  E(t). 

In  other  words,  the  clock  error  at  true  time  t  is  E(t),  and  all  other 
sources  of  error  in  TOA  measurement,  such  as  receiver  noise,  are  ignored 
here.  If  the  true  times  corresponding  to  Rq,  R^  and  R2  are  tQ,  t^  and  t2. 
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then  the  actual  transmission  time  that  would  correspond  to  reception  at  true 
time' t  will  be  given  by 


S 


o 


Sl  + 


(to  -  tl)  • 


if  range  acceleration  is  neglected.  It  is  easy  to  show  that  the  resulting 

A 

error  e  *  S  -  S  is  given  by 
o  o 


e  ■  E 

o 


where  E  *  E(t  ) „  etc,  and  the  factor 
o  o 


which  multiplies  the  entire  expression  for  e,  has  been  approximated  by  unity. 
This  is  permissible,  since  ds/dt  is  the  user-satellite  range  rate,  divided  by 
the  velocity  of  light,  and  is  hence  only  a  few  parts  in  10^.  By  the  same 
argument,  t ^  -  t^  is  essentially  equal  to  A,  and  we  write  e  in  the  final  form 


e  -  E  -  (l-X)E.  -  X  E, 
o  1  i 


where 
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The  true  receive  times  are  then  t,  ■  t  -  XA  and  t-  ■  t  +(1-X)A,  so  that 
the  error  is  completely  parametrized  by  the  quantity  X,  as  shown  In  Fig.  9-1. 
where  the  ranging  bursts  are  shown  as  simple  pulses.  This  error  In  trans¬ 
mission  time  will  be  used  instead  of  pseudorange  error,  since  the  reception 
time  will  be  common  to  all  satellites,  cancelling  out  of  the  position  cal¬ 
culation,  and  the  sign  is  of  no  importance. 


- mJ 
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FIG.  9-1  DEFINITION  OF  X 


Each  pseudorange  error  is  of  the  form  given  above,  each  with  its  own 
value  of  X  and  corresponding  receive  times,  while  all  share  the  error  term 
Eq.  If  the  clock  had  only  a  deterministic  error,  consisting  of  a  bias  and  a 
rate  error,  then  E(t)  would  have  the  form  a  +  bt,  and  it  is  easily  seen  that 
each  Individual  pseudorange  error  is  zero  in  this  case.  Physically,  this  is 
so  because  linear,  proportional  interpolation  is  Independent  of  the  origin 
and  scale  of  the  axis  against  which  the  data  is  plotted.  In  our  problem  this 
has  the  deeper  consequence  that  pseudorange  errors  are  insensitive  to  the 
low-frequency  components  of  random  clock  noise.  In  fact,  we  shall  find  that 
they  are  more  tolerant  of  flicker  noise  than  of  white  noise  in  the  clock 
output,  as  defined  below. 
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This  brings  up  the  question  of  clock  noise  models,  and  we  follow  the 
usual  practice  in  the  frequency  standards  field  by  representing  the  clock  as 
an  oscillator  whose  instantaneous  phase  at  time  t  is 

<|>(t)  -  27TfQ  /C  [1  +  y (s)  ]  ds  , 

where  f  is  the  nominal  frequency  of  the  oscillator.  The  quantity  y  is  the 
normalized  instantaneous  frequency  error: 

y(t)  .  (4>(t)/2ir)  -fQ 

fo 

A  clock  driven  by  this  oscillator  will  be  used  to  time  an  interval  by  measuring 
the  accumulated  phase  difference  over  the  interval,  and  dividing  by  2irfQ.  In 
terms  of  our  previous  notation,  the  error  in  timing  the  true  interval  (t^,  t2) 
will  be 

C2 

E(t_)  -  E(t1)  ■  /  y (s)  ds  . 

£1 

The  "frequency  noise",  y(t) ,  is  modelled  as  a  stationary  random  process  with 
zero  mean  (a  constant  mean  corresponds  to  a  constant  frequency  offset,  which 
has  no  effect  in  our  problem) . 

Two  models  are  used  here  for  the  y(t)  process.  In  the  first,  y(t)  is 
taken  to  be  "white  noise",  with  a  constant  power  spectral  density,  while  in 
the  second,  the  spectral  density  is  assumed  to  vary  inversely  with  frequency, 
corresponding  to  "flicker  noise".  The  effects  of  white  noise  are  easily 
analyzed  in  the  time  domain,  as  will  be  done  here,  but  flicker  noise  requires  a 
tedious  analysis  in  the  frequency  domain,  and  only  the  results  will  be  quoted. 


Clocks  are  usually  characterized  by  the  behavior  of  their  Allan 
variance*,  or  ( t) ,  a  dimensionless  quantity  which  measures  normalized  fre¬ 
quency  stability  over  an  interval  of  length  T,  as  a  function  of  t.  If 
the  frequency  noise  is  white,  with  single-sided  power  spectral  density  Gq, 
then  it  is  known  that  the  Allan  variance  is 
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varying  inversely  with  T.  Since  the  clock  frequency  noise  is  white,  the 
errors  in  timing  non -over lapping  times  are  uncorrelated,  and  the  error 
E(t')  -  E(t)  has  the  mean  squared  value 


G 

[  E(t ' )  -  E(t)]2  -  -f-  |t'  -  t|. 


(The  error  function  E(t)  is  therefore  a  Wiener  process).  The  expression 
allows  us  to  evaluate  the  mean  squared  values  of  pseudorange  errors,  and 
also  their  correlations,  in  terms  of  Gq  and  hence,  ultimately,  in  terms  of 
Allan  variance. 

Flicker  noise  has  spectral  density  G^f-1,  and  the  corresponding  Allan 
variance  is  independent  of  x 


a2  (x) 

a 


-  2  In  2  G 


*D.  W.  Allan,  "Statistics  of  Atomic  Frequency  Standards,"  Proc.  IEEE,  Vol. 
54,  No.  2,  February  1977,  pp.  221,  230. 
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Thi^s  connection  again  allows  us  to  express  our  results  directly  in  terms  of 
Allan  variance,  after  an  analysis  in  the  frequency  domain.  It  should  be  noted 
that  white  and  flicker  noise  models  are  in  general  agreement  with  the 
Allan  variance  behavior  of  quartz  crystal  clocks  in  the  range  of  time 
intervals  between  tens  of  milliseconds  and  seconds  which  are  relevant  to 
our  problem. 

Returning  to  our  expression  for  pseudorange  error,  we  note  that  it  can 
be  written  in  the  form 


e  -  (1-X)  (E  -  E.)  -  X  (E,  -  E  ) . 

oi  &  o 


Since  the  intervals  (t^,  tQ)  and  (tQ,  t,)  do  not  overlap,  we  easily  find 
the  mean  squared  value,  for  white  noise: 


G  7  G  A 

[(1-X)“  X  A  +  X2  (1-X) A]  -  X  (1-X). 


In  terms  of  the  Allan  variance  at  f'A,  we  can  put 


e2  -  a2  (A)  •  A2  •  X  (1-X)  . 
a 


Notice  that  the  error  vanishes  when  X  equals  zero  or  unity,  which  makes 
sense,  because  interpolation  is  unnecessary  when  a  pulse  arrives  exactly  at 
the  measurement  time.  A  similar  calculation  yields  the  mean  of  the  product 
of  two  pseudorange  errors  e  and  e': 


SI 


ee  -  a“  (A)  •  A2{XX'  (1-X')  +  (1-X)  ( 1 -X ’ )  X}  , 

where  X  is  the  smaller  of  the  two  parameters,  X  and  X'. 

With  these  formulas  we  can  evaluate  the  moment  matrix  for  the  clock- 
related  errors  for  any  number  of  pseudoranges,  corresponding  to  the  satellites 
currently  being  used  for  location.  This  moment  matrix,  in  turn,  fixes  the 
contribution  of  clock  errors  to  find  position  uncertainty.  Unfortunately, 
such  a  moment  matrix  depends  in  a  complicated  way  on  all  the  X-parameters 
which  describe  the  detailed  positions  of  the  received  pulses  relative  to 
the  chosen  measurement  time.  This  measurement  time,  however,  bears  no 
special  relation  to  the  pulse  arrival  times  and,  averaged  over  many  position 
updates,  using  different  sets  of  satellites,  the  X-values  will  appear  ran¬ 
domly  distributed  over  their  range,  from  zero  to  unity.  Since  we  only 
require  the  effect  of  clock  errors  on  mean  squared  position  error,  it  is 
reasonable  to  average  out  these  X-parameters,  yielding  a  much  simpler,  but 
representative,  pseudorange  error  moment  matrix.  By  integrating  over 
X  and  X',  from  zero  to  unity,  we  easily  find  the  mean  values 
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The  averaging  over  the  X-parameters  is  now  Included  in  the  significance  of 
the  overbar. 

According  to  our  previous  discussion,  <?a(A)  *  A  Is  the  rms  clock 
error  in  timing  an  interval  of  length  A,  and  hence  might  be  compared  with 
the  rms  error,  aR,  in  measuring  time  of  arrival  due  to  receiver  noise  and 
quantization  effects.  We  would  like  to  be  able  to  say  that  a  clock  error, 
a  (A)  •  A,  of  one  nanosecond  is  equivalent  to  a  TOA  error  of  some  other 
number  of  nanoseconds  in  its  effect  on  user  position  accuracy.  However, 

TOA  errors  are  independent  among  pseudorange  measurements,  corresponding  to 
the  values 


2 

e 


« 

ee 


TOA  errors. 


The  desired  comparison  is  made  possible  by  comparing  the  moment  matrices 
in  the  two  cases  for  pseudorange  differences,  rather  than  pseudoranges 
directly.  User  position  can  always  be  computed,  in  the  case  of  N  satellites, 
from  N-l  pseudorange  differences,  one  satellite  being  chosen  as  reference 
for  all  the  others.  Let  eQ  be  the  error  in  the  reference  pseudorange,  and 
let  e  and  e'  be  two  other  pseudorange  errors.  The  (N-l) -dimensional  moment 
matrix  of  differences  will  have  diagonal  elements  of  the  form 

(e  -  eQ)2  =  A, 


and  off-diagonal  elements  of  the  form 


f r 


(e  -  e  )  (e'  -  e  )  =  B. 

o  o 


We  can  work  these  out  in  terms  of  the  quantities  we  already  know: 

~~2  2 

A  -  e  +  e  ^  -  2ee 

o  o 

——  — ■  2 

B  ■  ee  -ee  -  e’e  +  e 
o  o  o 

For  TOA  errors,  the  cross-product  terms  are  zero,  and  all  pseudorange  errors 
2 

have  variance  crR  ,  hence 


A  - 

B  -  aR2  -  A/2 

(the  diagonal  elements  are  all  equal  to  A;  the  off-diagonal  elements  equal 
to  B).  For  clock  errors,  using  our  averaged  values,  g2  is  the  same  for 
all  satellites,  and  <  is  the  same  for  all  pairs,  which  means  that 


A  -  2  (e2  -  ee') 

B  -  -  ee  +  e2  -  A/2  . 


This  general  result  holds  for  all  noise  models;  for  white  noise  we  have 
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&  al  <A)  a2 


1  2 

vs  <«  4 


Since  these  moment  matrices  (for  pseudorange  differences)  have  the 
same  form,  except  for  a  proportionality  constant,  clock  errors  will  behave 
exactly  like  TOA  errors  in  their  effect  on  position  accuracy,  regardless 
of  the  number  of  satellites  used  for  location  on  the  geometry  of  the  con¬ 
stellation  in  use. 

The  proportionality  between  these  two  error  sources  is  expressed  by 
the  relation 


<°R>. 


a  (A)  *  A 

.  _§ _ _ 


■  0.258  0  (A)  •  A  (white  noise) 


which  gives  the  rms  TOA  error  equivalent  to  the  effects  of  a  clock  Allan 
2 

variance  o  (A) . 
a 

The  entire  analysis  can  be  repeated  for  flicker  noise,  which  results 
in  more  complicated  formulas  for  g2  and  ee-  in  terms  of  the  X-parameters. 
When  these  expressions  are  averaged  over  the  parameters,  we  obtain  the 
proportionality  relation 


(a0)  ■  0.159  a.  *  A  (flicker  noise) 

K  e  a 


where  a  ,  in  this  case,  is  independent  of  A. 
a 
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»  Taking  white  noise,  the  worst  of  the  two  models,  and  the  frame  length 
of  .5  seconds  for  A,  we  see  that  a  clock  whose  short-term  stability  (apart 
from  bias  and  frequency  offset)  is  characterized  by  an  Allan  variance  of  one 
part  in  107  (this  is  actually  the  square  root  of  the  Allan  variance) ,  for 
time  intervals  of  the  order  of  one  second,  contributes  position  errors 
equivalent  to  TOA  errors  of  rms  value  13  nanoseconds.  Since  these  errors 
are  combined  on  a  "root -sum-squared"  basis  with  other  errors  of  the  order 
30  nsec  or  more,  it  may  be  concluded  that  the  clock  stability  requirements 
of  the  TDMA  system  is  a  few  parts  in  107,  for  one-second  time  Intervals. 
Although  the  discussion  was  carried  out  in  the  context  of  a  pulsed  TDMA 
system,  the  results  also  apply  to  a  tone  ranging  system,  if  we  consider 
only  errors  in  timing  due  to  the  clock. 

The  PRFDMA  system  presents  a  similar,  but  more  difficult,  analysis 
problem.  In  steady  state,  each  tracker  is  smoothing  over  the  whole  past 
of  pulse  arrivals,  spaced  A  seconds  apart  (A  being  a  representative  pulse 
repetition  period).  The  chosen  measurement  time  will  be  XA  seconds  after 
the  latest  pulse  arrival  time,  where  X  varies  for  each  satellite,  being  in 
all  cases  between  zero  and  unity.  The  transmission  time  error,  e,  for  one 
satellite  is  then  a  weighted  sum  of  the  clock  errors  at  all  previous 
arrival  times.  As  before,  it  can  be  shown  that  e  vanishes  for  a  deter¬ 
ministic  clock  error  consisting  of  bias  and  rate  offset. 

For  white  noise  it  is  possible  to  evaluate  2  and  '  in  terms  of 

e  ee 

A  and  the  tracker  parameters,  and  the  time-phase  parameters  X  and  X'. 

Again,  X  and  X'  are  averaged  out  to  yield  an  effective  TOA  ranging  error 
of  the  form 


(a.)  -  K  a  (A)  •  A  , 

Re  a 
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where  K  Is  given  in  terms  of  the  tracker  parameters  by  the  expression 


„2  _  2a  -  (a-3)  (8/4)  +  (82/5)  , 

*■  3o  (4  -  2a  -  8) 


values  of  a  and  g  corresponding  to  a  tracker  "integration  time"  of  about 
10  pulses  are 


a  -  0.34545  ;  6  -  0.05455, 


which  results  in  K  *  0.438.  For  smoothing  over  5  pulses,  the  corresponding 
K  is  0.363.  These  K-values  are  larger  than  these  of  the  two  point  inter¬ 
polation,  mainly  because  transmission  time  is  being  extrapolated  outside 
the  measurement  interval,  but  the  difference  is  small.  In  our  application, 
the  pulse  Interval,  A,  is  only  a  few  milliseconds,  and  this  results  in 
much  less  demanding  clock  requirements  than  in  the  TDMA  case,  far  overcoming 
any  difference  in  K-factors  (flicker  noise  effects  have  not  been  evaluated 
for  the  alpha-beta  tracker,  but  should  result  in  smaller  K-factors  than 
white  noise) . 

As  an  example,  taking  A  ■  4  msec  and  K  *  0.5,  we  have 


aa(A) 


KA  (Ve 


or 


a  (.004  sec) 

A 


6.5  x  10'6  . 


To  be  fair,  this  same  clock  will  have  a  smaller  Allan  variance  at  T  *  .5  sec 
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(the  specification  point  for  the  TDMA  case  discussed  above),  with  our  white 
noise  model,  namely 


o  (.5 

a 


sec) 


"(soo) 


x  6.5  x  10 


-6 


-  5.8  x  10 


-7. 


hence,  the  clock  requirements  are  relaxed  here  only  by  a  factor  of  about  five. 


APPENDIX  A 


ANALYSIS  OF  PERFORMANCE  OF  GPS  C/A  WAVEFORM  WITH  MORE  POWER 

This  appendix  provides  some  analytical  details  to  back  up  the  summary 
statements  of  Section  3.  The  outline  of  this  appendix  is  the  same;  first  we 
discuss  code  synchronization  issues,  then  carrier  phase  lock  issues  and  finally 
data  demodulation  and  issues  related  to  data  gathering. 

A.l  Code  Synchronization 

A. 1.1  Sequential  Code  Synchronization 

The  average  time  required  for  code  synchronization  is  a  function  of  the 
assumed  signal- to-noise  power  density.  This  study  assumes  an  RF  transmitted 
power  by  the  GPS  satellites  of  from  14  dBW  (present  level)  to  29  dBW.  Table  A-l 
summarizes  our  estimates  of  link  losses  and  the  ratio  of  signal  power  to  noise 
power  density  at  the  input  of  the  receiver. 

The  two  issues  to  be  considered  are:  the  time  required  for  sync  detec¬ 
tion,  and  the  sensitivity  to  doppler.  If  the  system  is  sensitive  to  doppler 
then  one  must  not  only  perform  a  code  search,  but  a  search  in  the  frequency 
dimension  as  well.  In  a  serial  code  search,  the  replica  code  is  tested  against 
the  incoming  code  with  steps  of  one  half  chip  length.  Detection  is  based  on 
the  power  level  of  the  despread  carrier  falling  in  the  narrow  bandpass  filter 
of  the  sync  detector.  With  Tg  denoting  the  signal  integration  time,  it  is 
assumed  that  the  incoming  signal  is  not  doppler  shifted  in  frequency  as  com¬ 
pared  to  the  local  oscillator*  by  more  than  2Y  This  frequency  offset 

causes  a  loss  of  no  more  than  3.8  dB  in  output  signal  to  noise  ratio  at  the 
correlation-peak.  Another  1.25  dB  is  lost  in  case  the  two  codes  are  one 
quarter  chip  misaligned  and  the  testing  is  done  with  one  half  chip  increments. 

The  slgnal-to-noise  ratio  must  be  15  dB,  i.e., 

C  T  >15  dB  , 

N  3 
o 

to  guarantee  10  dB  under  the  worst  conditions  of  maximum  frequency  shift  and 
code  misalignment.  The  10  dB  corresponds  to  a  probability  of  false  alarm  of 
less  than  one  percent  combined  with  a  probability  of  correct  detection  of  90Z. 
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TABLE  A-l.  HIGH  POWER  GPS  BUDGET 


RF  Sat.  Transmitted  Power 

Sat.  Antenna  Gain 

Pathless 

Atmospheric  Loss 

14  to  29  dBW 

11.5  dB 

-185.5  dB 

-1  dB 

At  Receiver  Antenna 

Receiver  Antenna  Gain 

Receiver  Loss 

Margin  for  Other  Loss 

-161  to  -146  dBW 

0  dB 

-2.5  dB 

-2  dB 

-165.5  to  -150.5  dBW 

Receiver  Noise  for  5  dB  Noise  Figure 

-199  dBVtfHz 

C/N  ;  Ratio  of  Carrier  Power  to 
o 

Noise  Power  Density 

33.5  to  48.5  dBHz 

r 

Since  ^  is  variable  from  33.5  Co  48.5  dBHz,  Che  incegration  time  ranges  from 
14. 3„  to0. 45  msec.  Since  the  code  epoch  of  the  presently  proposed  GPS  C/A 
signal  is  1  msec  long,  one  will  prefer  coherent  integration  times  which  are 
multiples  of  1  msec.  If  a  different  GPS  signal  format  with  shorter  epochs 
were  used  in  combination  with  more  power,  coherent  integration  times  less  than 
1  msec  could  be  exploited.  We  expand  on  this  in  a  later  paragraph. 

The  above  integration  times  obviously  depend  on  the  threshold  set  for 
reliable  sync  detection.  If  a  higher  degree  of  false  alarm  and  missed  detec¬ 
tions  is  acceptable,  shorter  times  will  result. 

It  would  generally  be  unwise  to  have  coherent  integration  times  of  over 
19  mseconds  to  reach  the  required  detection  energy,  as  would  be  necessary  at 
the  lower  scale  of  the  power  level  since  this  would  reduce  the  filter  band¬ 
width  to  less  than  75  Hz.  This  makes  the  code  synchronization  function  much 
more  sensitive  to  doppler  and  frequency  drift  from  the  oscillator. 

It  would  not  be  useful  to  have  integration  times  less  than  1  msec  for  the 
present  C/A  code  since  appropriate  threshold  levels  for  partial  code  auto 
correlations  are  difficult  if  not  impossible  to  set  reliably. 

A  coherent  integration  time  of  1  msec  allows  roughly  a  frequency  band  of 
±  500  Hz.  For  a  sequential  receiver  cycling  through  four  satellite  signals 
with  cycle  times  on  the  order  of  1  to  10  seconds,  the  frequency  drift  accumu¬ 
lated  from  relative  user  and  satellite  motion  and  clock  drift  will  generally 
be  well  within  the  ±  500  Hz  band.  This  avoids  a  sequential  code  sync  search 
over  several  frequency  bands. 

The  situation  is  entirely  different  when  the  receiver  is  trying  to  acquire 
the  signal  from  a  new  satellite.  The  number  of  frequency  bands  to  be  searched 
depends  on  the  a  priori  frequency  uncertainty.  Even  with  a  perfect  oscillator, 
the  carrier  frequency  can  be  off  as  much  as  ±  4342  Hz  (worst  case  at  the 
earth  poles)  from  satellite  motion  with  respect  to  a  stationary  user  and  ±  1000 
Hz  from  user  motion.  If  sufficiently  accurate  satellite  ephemeris  data  are 

available,  for  example,  in  the  form  of  pre-stored  almanac  data,  only  the  latter 

\ 

\ 
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frequency  error  is  important,  implying  a  search  in  at  most  2  or  3  frequency 
steps.  If  no  almanac  data  are  available,  as  would  be  the  case  in  a  "cold 
start-up”  mode,  as  many  as  10  frequency  steps  may  be  necessary. 

Since  we  have  assumed  that  the  code  synchronization  search  is  done  sequen¬ 
tially  in  one-half  chip  steps,  it  could  require  up  to  2046  steps.  The  inte¬ 
gration  time  discussed  above  must  then  be  multiplied  by  this  number  and  by  the 
number  of  required  frequency  steps  which  can  be  from  2  to  10,  depending  on 
power  level. 

In  summary,  the  effect  of  more  power  on  sequential  code  synchronization 
is  that  integration  time  (inversely  proportional  to  power)  can  be  significantly 
reduced,  but  has  a  practical  lower  limit  of  one  msec  (code  epoch).  The  number 
of  frequency  steps  can  also  be  reduced,  since  that  number  equals  the  assumed 
frequency  uncertainty  present,  divided  by  the  bandwidth  (inverse  of  integration 
time) .  While  the  reduction  can  be  dramatic  for  start-up  conditions  when  the 
code  is  acquired  the  first  time,  the  number  of  frequency  steps  is  more  likely 
to  be  two  in  normal  operations  with  no  reductions  due  to  increased  power 
possible  when  cycling  times  between  satellite  signals  are  from  1  to  10  seconds. 
The  largest  factor  in  synchronization  is  the  number  of  code  phases  to  be 
tested.  That  number  is  not  directly  dependent  on  power  level,  certainly  not 
in  the  acquisition  phase.  However,  in  the  re-acquisition  mode  (when  cycling 
between  satellite  signals)  the  number  of  code  phases  to  be  searched  is  a 
function  of  the  range-uncertainty  accumulated  during  the  time  other  satellites 
were  visited.  More  power  can  provide  better  ranging  (see  Section  3.1.4)  and 
more  accurate  range  prediction  and  therefore  will  reduce  the  number  of  code 
phases  to  be  searched  in  the  re-acquisition  mode. 

The  above  analysis  points  to  shorter  synchronization  times  with  increased 
power.  This  is  an  important  operational  benefit.  Hidden  benefits  are  in  the 
potential  for  use  of  a  cheaper  IF  oscillator,  elimination  of  software  to 
perform  incoherent  integration  and  eventually  elimination  of  software  associated 
with  frequency  searching  in  normal  operations.  Frequency  searching  would  be 
done  only  when  a  new  satellite  is  being  acquired. 

The  possibility  of  elimination  of  code  phase  searching  is  suggested  in 
Section  A.  1.3  on  matched  filtering. 


A. 1.2  Possible  Use  of  a  Shorter  Code 


the  previous  section  suggests  that  if  available  signal-to-noise  power 
density  ratio  exceeds  45  dBHz,  then  with  the  same  detection  thresholds  as  above 
the  integration  time  for  the  test  on  one  code  phase  in  one  frequency  band 
could  be  less  than  1  msec.  The  objection  against  lesser  time  was  that  the 
code  epoch  was  1  msec  and  constituted  a  lower  bound.  That  objection  is  removed 
if  one  is  willing  to  redesign  the  GPS  waveform  with  a  shorter  code  of,  let's  say, 
511  or  even  127  chips. 

The  potential  of  the  short  code  of  127  chips  at  the  present  1.023  Mbps 
clocking  rate  is  that  one  can  coherently  integrate  over  one  epoch  and  add  the 
results  incoherently  to  reach  a  sufficient  energy  level  for  detection  and  be 
totally  independent  of  the  maximum  expected  doppler.  Tolerable  frequency  off¬ 
sets  are  ±  7900  Hz  for  a  127  chip  code.  Circuitry  and  algorithms  corresponding 
to  frequency  stepping  can  be  dropped  and  a  low  quality  oscillator  is  permissible. 

But  all  this  has  a  price.  Shorter  codes  result  in  a  lower  margin  between 
the  auto  correlation  peak  and  the  cross  correlation  peaks.  That  margin  is 
difficult  to  calculate  except  for  Gold  codes.  The  margin  can  also  be  affected 
by  the  carrier  frequency  offset  existing  between  the  cross  correlated  signals 
if  that  offset  is  larger  than  several  inverse  integration  times.  For  example, 
for  a  1023  bit  long  Gold  code  with  1  msec  integration  time,  doppler  offsets  of 
several  kHz  among  the  signals  of  the  (at  most)  eleven  satellites  do  exist. 

The  margin  discussed  above  is  in  the  worst  case  about  3  dB  lower  than  the  mar¬ 
gin  given  by  the  Gold  bound  expressing  maximum  cross  correlation  power  P  at 
zdro  frequency  offset: 


Obviously,  given  a  set  of  24  Gold  codes,  the  exact  margins  for  the  appli¬ 
cable  frequency  offset  range  can  be  obtained. 


For  a  code  as  short  as  127  chips  and  corresponding  short  integration 
times,  no  applicable  doppler  offset  can  exceed  several  inverse  times  and 
the  above  effect  does  not  exist.  The  margin  expressed  by  the  Gold  bound 
then  holds  for  all  pairs  of  satellite  signals;  for  a  seven-register  code, 
P  equals: 


On  the  other  hand,  insensitivity  to  doppler  means  that  now  more  satellite 
signals  are  capable  of  interfering.  Since  as  many  as  eleven  can  be  simultan¬ 
eously  visible  more  signal  interference  exists.  Several  dB  in  margin  can  be 
lost  due  to  multiple  signal  interference. 

In  summary,  the  negative  aspects  of  using  a  shorter  code  are  that  the  mar¬ 
gin  between  auto  correlation  peak  and  cross  correlation  peak  is  lowered.  But, 
even  for  a  short  127  chip  code,  a  healthy  17.5  dB  margin  continues  to  exist. 

Even  with  more  signals  being  capable  of  interfering,  multiple  code  accessing 
capability  is  sufficient. 

We  conclude  this  section  on  shorter  codes  with  some  comparison  curves  for 
code  acquisition  in  the  worst  case,  i.e.  all  code  phases  must  be  searched  and  all 
frequency  bins  must  be  searched  for  the  maximum  doppler  that  possibly  can 
occur  and  the  threshold  for  detection  of  synchronization  is  a  guaranteed  10  dB 
(even  with  the  code  time  misalignment  of  one  quarter  chip  and  the  frequency 
offset  between  incoming  signal  and  locally  generated  signal  equal  to  one-half  the 
frequency  band  or  equal  to  the  inverse  of  the  coherent  integration  time) . 

Figure  A.-1  presents  the  time  required  for  a  detection  decision  on  a  single 
code  phase  for: 

a)  Curve  A:  A  short  127  chip  code  and  incoherent  sequential  de¬ 
cision.  The  coherent  integration  time  is  taken  to  be  about  125  ysec  (corre¬ 
sponding  to  a  127  chip  code  length  at  1.023  MHz  clock  rate).  The  result  of 
each  integration  is  the  sum  of  in-phase  and  quadrature  component  energy  and 
is  added  N  times.  The  number  N  is  chosen  to  represent  a  final  probability  of 
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_3 

detection  of  0.9  and  a  probability  of  false  alarm  of  5  x  10  .  The  average 
value  for  N,  denoted  N,  can  be  derived  from  radar  considerations  and  is  given 
roughly  by 

£  ® 

where  ^  is  the  ratio  of  energy  per  pulse  from  coherent  integration  to  the  one- 
sided  nBise  spectral  density  expressed  in  dB.  For  example,  for  C/N  -  43.8  dB 
Hz  (which  is  10  dB  above  the  present  GPS  signal  level)  and  Tj  -  0.125  msec, 

E/Nq  ■  4.8  dB  and  N  ■  6.  Effectively  one  needs  at  least  one  more  step  to  ex¬ 
ceed  the  threshold.  Plotted  is  therefore  Tj  (N  +  1)  ■  0.9  msec. 

b)  Curve  B  corresponds  to  the  coherent  integration  time  needed 
to  obtain  the  required  E/Nq  when  a  filter  bank  is  used  with  as  many  frequency 
bins  as  needed  in  parallel  to  cover  the  total  doppler  D  assumed  to  be  ±  4  kHz 
(user  motion  induced  doppler  and  maximum  doppler  at  mid  latitude) .  The  required 
E/Nq  was  held  at  13  dB  to  obtain  similar  detection  and  false  alarm  capabilities. 

c)  Curve  C  is  the  total  test  time  if  only  one  filter  exists  which 

sequentially  tests  over  all  frequency  bins.  If  Tj*  is  the  coherent  integration 
time  the  total  time  Tg  equals 


d)  Curve  D  simply  represents  the  number  of  frequency  steps 

corresponding  to  curve  B  and  is  inversely  proportional  to  assumed  C/N  . 

o 

^  clear  that  short  incoherent  sequential  detection  is  comparable  to 
detection  with  coherent  multiple  velocity  filter  detection  as  far  as  detection 
time  on  a  single  code  phase.  The  sequential  nature  of  the  detection  is  traded 
against  multiplicity  of  filters  albeit  with  a  slight  loss  in  efficiency. 

Figure  A-2  presents  again  curve  A  of  Fig.  A— 1  compared  to  curve  E  repre¬ 
senting  the  more  realistic  GPS  case  of  coherent  integration  of  over  1  msec  and 
incoherent  sequential  integration  to  satisfy  the  E/Nq  requirement  of  13  dB 
combined  with  the  use  of  a  single  velocity  filter. 
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TEST  TIME  (nsec) 


FIG.  A-l  TEST  TIME  ON  ONE  CODE  PHASE  FOR  INITIAL  CODE  ACQUISITION  FOR 
MAXIMUM  FREQUENCY  UNCERTAINTY  OF  +  4  kHz 
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NUMBER  OF  PARALLEL  FILTERS  FOR  CURVE  B 


Test  Time  (msec) 


FIG.  A- 2  TEST  TIME  ON  ONE  CODE  PHASE  FOR  INITIAL  CODE 
ACQUISITION  FOR  MAXIMUM  FREQUENCY  UNCERTAINTY 


It  is  clear  that  short  codes  are  advantageous  for  code  acquisition  and 
for  simplicity  of  acquisition  circuitry  due  to  their  independence  of  doppler. 

The  cost  is  a  possible  false  lock  due  to  the  reduced  cross  correlation  margin 
and  higher  sidelobes  of  the  auto  correlation  function.  Also,  the  range  ambi¬ 
guity  Interval  is  smaller  for  shorter  codes  and  may  require  more  software 
manipulations . 

A. 1.3  Use  of  Matched  Filter 

If  instead  of  the  sequential  digital  filter  one  uses  a  matched  filter  such 
as  a  CCD  correlator,  one  does  not  have  to  test  every  code  phase  in  one  half  chip 
increments  and  there  is  a  corresponding  saving  in  acquisition  time.  Programmable 
CCD  devices  are  available  for  code  lengths  of  127.  Programmable  CCD  devices 
for  code  lengths  as  long  as  1023  will  not  be  available  in  the  near  future. 

The  short  code  and  sequential  incoherent  detection  method  therefore  can  benefit 
from  this  technology  to  great  advantage. 

The  whole  operation  not  only  of  code  synchronization  but  also  of  code 
tracking  the  PLL  and  demodulation  could  be  done  at  baseband. 

A. 1.4  Improved  Code  Tracking 

Once  synchronization  is  established ,  despreading  of  the  signal  ce  i  be  done 
for  the  purposes  of  demodulation  of  the  data  and  obtaining  a  doppler  reading. 

A  feedback  from  the  PLL  "tunes"  the  locally  generated  signal  to  maximize  the 
error  readings  for  any  choice  of  discrimination  curve.  These  error  readings 
are  also  a  function  of  signal-to-noise.  The  chip  duration  A  can  be  "split" 
more  accurately  with  Increased  C/Nq  or  in  other  words  more  accurate  range 
measurements  can  be  expected. 

Ranging  error  expressed  as  the  variance  of  chiplength  A  in  terms  of  C/Nq 
is  (Ref.  9). 


where  k  ■  0  for  coherent  tracking;  ■  1  for  a  non-coherent  DLL 

is  one-sided  bandwidth  of  the  code- tracking  loop 
Bv  is  the  inverse  of  the  Integration  time  of  the  predetection 
bandpass  filter  in  the  non-coherent  DLL. 
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This  equation  shows  for  example,  that  a  10  dB  increase  in  C/Nq  gives  a  chip 
splitting  capability  at  least  three  times  better. 


A.  2  Performance  of  Carrier  Tracking  Loop 


The  PLL  is  used  principally  to  demodulate  the  data  coherently.  A  higher 
C/Nq  level  provides  better  performance  and  allows  larger  bandwidth.  The  latter 
makes  the  system  less  sensitive  to  relative  user /satellite  dynamics.  For  a  2nd 
order  loop  the  acceleration  lag  6^  1x1  response  to  a  constant  acceleration  step 

A&  (rad/sec^)  is 


9  m  Aw 

Lag  2 

°  w 


For  a  critically  damped  loop  the  loop  bandwidth  relates  to  the  eigenfrequency 

w  by  1  ■  0.53  w  .  The  noise  term  contribution 
in  L  in 


It  can  be  seen  that  the  value  E  of  total  phase  error 


k0T  + 
Lag 


will  be  optimized  for  different  based  on  the  weighting  k  and  the  assumed 

dynamics,  but  is  strongly  affected  by  the  assumed  C/N  .  Increased  C/N 

o  o 

reduces  the  phase  error  E  and  allows  the  loop  bandwidth,  which  is  optimized 
for  the  noise  contribution  and  the  lag  component,  to  be  larger.  This  is  shown 
more  clearly  in  Fig.  3  1  of  Section  3  for  the  indicated  weighting.  The  choice 
of  weighting  is  obviously  up  to  the  designer.  An  acceptable  choice  for  k  is  2. 
allowing  the  PLL  to  cope  well  with  dynamics  as  high  as  several  g. 

A. 3  Data  Demodulation 


High  C/Nq  will  allow  better  data  demodulation.  First,  the  larger  bandwidth 
for  the  Costas  loop  will  mean  fewer  phase  slips  of  that  loop  when  subjected  to  a 
dynamic  situation  while  maintaining  a  small  phase  tracking  error.  The  cosine  of 
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Chat  error  affecCs  Che  detection  energy  per  bit  collected  over  20  msec.  The 
energy  is  obviously  proportional  to  C/Nq  which  can  then  be  directly  translated, 
for  a  DPSK  data  modulation,  to  bit  error  rates. 

A. 4  Satellite  Scheduling 

In  the  previous  sections  the  sensitivity  of  acquisition  and  re-acquisition 
time  to  SNR  was  demonstrated.  In  a  sequential  receiver  the  latter  depends 
also  on  cycle  time.  Let  C  denote  the  cycle  time  during  which  the  four  satellites 
selected  for  user  position  determination  have  been  "visited". 

For  shorter  C,  less  range  uncertainty  to  a  particular  satellite  will  have 

accumulated  while  the  signals  from  the  three  others  are  being  tracked.  In 
the  software  algorithms  range  and  range  rates  are  assumed  to  be  tracked. 

Obviously  shorter  C  means  higher  track  update  rate  and  correspondingly  more 
accurate  estimates  of  ranges,  range  rates,  clock  bias,  user  position  and 
velocities  due  to  the  smoothing  process  applied  to  the  data.  The  uncertainties 
on  predicted  range  (predicted  code  delay)  and  range  rate  (doppler  offset)  are 
then  mostly  determined  by  unexpected  accelerations. 

Accelerations  from  satellite  motion  are  less  than  one  wavelength  per 
second  squared.  The  biggest  accelerations  derive  from  user  motion  and  are 
limited  to  a  few  tens  of  meters  per  second  squared.  Range  uncertainty  accumulated 
over  the  three-quarter  cycle  time  C  therefore  will  seldom  exceed  one  chip 
length  (300m)  for  C  up  to  5  seconds. 

Therefore  re-acquisition  may  require  only  the  checking  of  one  to  three  half  chip 
code  phases.  Fig.  A-l  indicated  the  test  time  for  high  SNR  still  only  to  be 
of  the  order  of  tens  of  msec.  Accelerations  also  cause  frequency  offset  to 
accumulate  linearly  with  cycle  time  during  the  three-quarter  cycle  time.  To 
this  one  must  add  the  oscillator  error  during  the  three-quarter  cycle  time. 

For  a  master  oscillator  of  10  ^  Hz /Hz  the  frequency  error  accumulates  at  the 
rate  of  160  Hz/sec.  The  DLL  for  code  synchronization,  for  high  SNR,  can  be 
designed  with  a  wide  enough  bandwidth  (e.g.,  1  kHz)  so  that  no  frequency  steps 
must  be  considered.  The  PLL  will  then  require  AFC  aiding  to  cope  with  such 
large  initial  frequency  offsets. 


What  is  saved  on  oscillator  cost  therefore  may  well  have  to  be  paid  in  a 
more^ complicated  PLL  unless  one  reduces  the  cycle  time  drastically  to,  for 
example,  two  seconds.  The  disadvantage  here  is  that  the  effective  time  for 
acquiring  data  becomes  very  small  and  message  assembly  of  out-of-sequence 
portions,  with  no  guarantee  of  message  completeness  due  to  missing  data  bits, 
becomes  difficult. 

Since  tracking  and  data  collection  must  occur  simultaneously,  the  better 
approach  is  for  longer  cycle  times,  and  a  better  oscillator  (e.g.,  10  Hz/Hz). 

Another  point  to  be  emphasized  is,  that  if  ranging  accuracy  can  be  relaxed 
the  DLL  which  provides  the  ranging  potential  can  be  simplified.  Instead  of 
using  parallel  channels  to  provide  early/late/on  time  outputs  of  the  local 
code  generator,  a  well  known  timesharing  method  called  tau  dither  can  be  used 
with  only  slight  degradation  in  performance. 


APPENDIX  B 


DETAILS  OF  TDMA  PULSED  GPS  RECEIVER 

In  Section  4'  a  design  was  described  which  represents  the  class  of  TDMA 
Pulsed  waveforms.  Several  small  studies  were  performed  to  determine  the 
feasibility  and  simplicity  of  the  approach. 

In  this  Appendix  we  include  some  details  on  the  sync  processing  and  data 
processing  logic  and  how  the  logic  translates  into  hardware.  The  flow  charts 
and  diagrams  are  otherwise  self  explanatory. 

Figure  B-l  presents  the  processing  waveform.  Figure  B-2  shows  the  sync 
processing  flow  chart  translating  into  logic  cards  D^  to  Dj  in  Figures  B-3  to 

B-5.  Figures  B-6  shows  the  data  demodulator  counter  states  and  Figure  B-7 
summarizes  the  layout  of  the  GPS  receiver. 
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SPECIAL  PURPOSE  LOGIC  CARD  D-l 


FIG.  B-4  SPECIAL  PURPOSE  LOGIC  CARD  D-2 


FIG.  B-5  SPECIAL  PURPOSE  LOGIC  CARD  D-3 
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FIG.  B-7  PULSED  GPS  RECEIVER:  CASE  REMOVED 


APPENDIX  C 


POWER  DISTRIBUTION  IN  TONES  FOR  GPS  TONE  RANGING  SYSTEM 

In  che  cone  ranging  syscem  described  in  Seccion  5  Che  FM  or  PM  modulation 
index  is  chosen  so  chac  Che  power  in  che  fundamencal  speccral  lines  is  maxi¬ 
mized.  In  che  cone  ranging  mechod,  no  convenCional  FM  or  PM  demodulacion  is 
used.  Inscead,  a  PLL  is  applied  Co  recover  Che  specific  Cone,  after  heterody- 
ing  che  modulaced  carrier  Co  baseband.  ThaC  allows  us  Co  realize  an  increase 
of  3  dB  in  loop  SNR  over  che  SNR  ac  IF.  Ic  is  imporcant  Co  mainCain  sufficienc 
power  in  che  carrier  frequency  speccral  line  during  Cone  modulacion  in  order  to 
allow  the  carrier  PLL  to  track  with  sufficiently  small  phase-error.  The  signal 
amplicude  of  Che  cone  is  proportional  to  che  cosine  of  that  phase  error. 

In  this  appendix  we  derive  an  expression  for  the  SNR  in  the  tone  PLL  as  a 

function  of  modulation  index.  The  index  is  then  chosen  to  maximize  the  ratio 
T 

77  ,  signal- to-noise  power  density  ratio  at  the  input  of  the  tone  PLL  to  the 
C 

ratio  —  ,  signal- Co-noise  power  density  at  the  input  of  the  receiver, 
o 

In  the  following  derivation  C(t)  will  represent  the  received  RF  waveform, 
the  set  f'  are  the  tones,  f£  is  the  carrier  frequency,  n(t)  is  the  noise  com¬ 
ponent  of  Che  RF  signal. 

Let  the  carrier,  frequency  f£  ,  be  phase  modulated  by  the  sum  of  N  tones 
f '  .  The  received  RF  waveform  is  then 

C(t)  -  A  cos  f 2irf  t  +  <f>  +  E  A  •  sin(2irf'  •  t)l  +  n(t) 

where  n(t)  -  nc(t)cos (2irfct)-ng(t)sin ^irf^t) 

The  noise  components  n  and  n  are  low  pass,  gaussian  with  zero  mean,  and 

c  s 

standard  deviation  equal  Co  standard  deviation  of  n.  If  n  has  single-sided 

spectral  density  N  in  che  RF  band  B,  Chen  n  and  n  have  density  2N  in  lowpass 
g  o  cs  o 

band  .  (Nq  •  kTN^,  ,  Nj,  is  noisef  igure) .  In  terms  of  the  received  waveform 
the  signal- to-noise  power  density  ratio  ^  equals: 
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C  m  A 
N  “  2N 


N 

Q  A  A  7T 
i-1 


<A±) ,  where  ^  are  Bessel  functions  of  i 


th 


order , 


N  J,(A.) 

then  C(t)-n(t)  -  Qcos(2irf  t  +  $)  +  Q  E  -7~— A-  cos  [2ir(f  +  £'  )t  +  $] 

C  i-1  W  c  1 


4  J.(A  ) 

-  Q  Z  T~7T  V  cos  [2irf  -fpt  +  $1  + - 

i-1  W  C  1 

higher  order  terms. 

The  higher  order  terms  are  harmonics  and  intermodulation  products  at 
spectral  frequencies 

N 

fc  v,  v; 

i-1 

for  all  possible  combinations  [k^  k2>  k-j,...^]  for  positive  and  negative 
integers  k^ 

Assuming  that  an  accurate  carrier  reference  is  obtained  during  acquisition, 
a  quadrature  reference 


b(t)  -  a  sin  (2 it  f  +  <J>) 

A  /V 

can  be  formed.  For  a  PLL  fc  -  fc,  tp  -  <f>  +  5<f>,  6$  being  the  carrier  phase 
tracking  error. 

The  received  signal  C(t)  is  heterodyned  to  baseband  by  this  quadrature 
reference  b(t).  Let  C'(t)  -  C(t)-n(t);  then 


<  C'(t)  b(t) 


> 


Low  Pass 


w 

w 


sin  (2irf '  t-<S4>) 


111 


N  J.  (4, ) 

+  f9,J,  T7S7T  3ln(-2,,fI  t"s*) 

£  1"X  O  1 


N  J,(A,) 

*  -  aQ  Z  ~  yr-v  .  cos  64*  sin  (2irfjt). 

1*1  o^ai/  1 

Observe  that  Che  amplitude  of  the  Cone  fundamental  at  baseband  Is,  for  tone  f'  , 

vv 

cos  6<j>.  aQ  J—  «  aJCAj^)  .JQ.  (A2)  . .  .  cos  54> 

o '  1 

We  now  discuss  Che  term  n(t)  heterodyned  to  baseband  : 

<  n(t)  b(t)>Low  pggg  -  f  (ncsin<{)  -ng  cos  0) 

B  2  2 

This  is  low  pass  noise  over  bandwidth  ■=■  and  spectral  density  ~  2  N  ■  a  --■a  - 

£  4  o  2 

The  tone  f*  is  now  isolated  in  a  PLL  with  tone-to-spectral  noise  density ( 

of 


I  (aQ  J1  (V\cos2gO 

(T  \  4  2  \  Jo  (V  /  A2  J2(4. )..J2<4 

W  2  „  50  »  1  1  " 

n  a  No 


Jq  (An)  cos250  , 


.  C  A 
or,  since  -  Jn 

o  o 


(i)-2f  Jo  (41)..jJ(4n)...J2(V=os2S* 

'  °'n  o 


A  similar  analysis  shows  that  the 


N 

C  v  j  (A.)»cos  60  . 

N  i-1  0  1 

o 


(0 


o'  carrier 


at  the  RF  frequency  f  is 
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This  is  the  power  ratio  available  to  the  carrier  PLL  for  carrier  tracking  and 
oust  be  sufficient  to  maintain  accurate  phase  track  since  the  cosine  of  the 
phase  error  scales  the  power  available  to  the  tone  PLL. 

The  method  of  clustered  tones  allows  the  range  accuracy  to  be  determined 
from  the  average  of  all  tone  phases  while  ambiguity  resolution  is  based  on  the 
sum  or  difference  of  two  tones  implying  a  range  error  which  is  a  factor  /3T~ 
worse  chan  the  error  associated  with  a  single  tone.  This  effectively  favors 
the  finest  of  the  fundamental  set  of  tones  compared  with  all  other  tones  and 
is  a  desirable  feature.  For  simplicity  one  may  allow  the  clustered  tones  to 
have  equal  modulation  indices  A. 

For  any  tone  PLL  one  may  then  write 


where 

2  2N'2  2 

k  -  2  cos  Jq(A)  (A)  . 

For  6<fr  «  o  and  N  ■  4,  Fig.  5“  3  in  Section  5  shows  k  as  a  function  of  modulation 
index  A.  Factor  k  is  sometimes  called  the  suppression  factor.  The  optimum  A  is 
.7  for  minimum  suppression.  However,  as  can  be  seen,  the  minimum  is  fairly  flat 
for  A  between  0.5  and  0.9.  One  may  favor  smaller  A  to  limit  the  maximum  phase 
deviation  (2.8  radians  for  4  tones  and  A  ■  0.7)  or  limit  the  effective  band¬ 
width  which  is  occupied.  The  approximate  bandwidth  of  the  RF  signal  is  given 
by  the  Carson  rule  , 

2(A1  +  l)f1  , 

since  all  tones  are  clustered  around  and  below  f. .  The  SNR  in  the  tone  PLL, 

c  ^ 

for  rr  equal  to  the  projected  45  dBHz,  and  A  *  .7  would  equal  38.1  dBHz.  The 
wo 

design  assumed  35  dBHz  to  allow  about  3  dB  of  margin. 
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